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Preface 

Vheu  the  Air  Cushion  System  replaced  the  conventional  takeoff 
and  landing  systems  of  the  Jlndlvik  remotely  piloted  vehicle,  the 
possibility  existed  that  Instabilities  in  pitch,  roll,  and  yaw 
could  occur.  As  a result,  this  paper  was  intended  as  a design 
of  a Takeoff  Stabilization  System  for  the  Jlndlvik  using  existing 
autopilot  sensors  and  Incorporating  an  engine  yaw  thruster  and 
vertical  wing  tip  roll  thrusters.  When  the  design  was  completed, 
it  was  sufficiently  general  that  the  technique  could  be  applied 
to  any  air  cushion  aircraft  or  VTOL  aircraft.  The  Landing  Sta- 
bilization System  for  the  Jlndlvik  using  the  same  sensors  and 
actuators  Is  presently  being  designed  by  Captain  Max  Stafford 
as  his  thesis  for  the  Air  Force  Institute  of  Technology  (AFIT) . 

' 1 wish  to  express  my  gratitude  to  my  Thesis  advisors.  Dr. 

George  Kurylowich  of  the  Air  Force  Flight  Dynamics  Laboratory 
(AFFDL)  and  Major  R.  Potter  of  AFIT.  Also,  thanks  are  due  to 
Captain  James  Negro  of  AFIT,  Major  Jack  Randall  and  Mr.  Jim 
Steiger  of  the  Air  Force  Flight  Dynamics  Laboratory  for  their 
technical  advice  and  assistance. 

My  wife,  Jane,  does  not  knew  how  much  she  has  contributed 
to  this  study,  but  her  patience,  understanding,  and  encouragement 
has  definitely  made  the  past  eighteen  months  of  work  much  easier. 

Edward  Kenney 
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Abstract 

The  Inherent  Instability  In  pitch  and  roll  associated  with 
an  Ai.r  Cushion  Landing  System  (aCLS)  aircraft  at  low  airspeeds 
was  investigated,  and  a means  to  aid  control  in  pitch  and  roll 
was  developed.  The  control  system  required  the  use  of  vertical 
wing  tip  thrusters  which  provided  thrust  up  or  down  depending  on 

the  control  signal  (similar  to  space  vehicle  thrusters) . These 

« 

thrusters  could  be  activated  alternately  to  control  roll  angle 
and  roll  rate  with  the  use  of  a bang-bang  optimal  controller. 

As  well,  the  thrusters  would  be  set  forward  of  the  aircraft 
centre  of  gravity  and  could  be  activated  in  tandum  to  aid  in 
pitch  control. 

The  Jindivlk  Remotely  Piloted  Vehicle,  an  Australian  tar- 
get drone,  was  fitted  with  an  ACLS  and  taxi  tests  showed  the 
Instability  and  need  for  a stabilization  system.  Subsequent 
use  of  Jindivlk  wind  tunnel  and  taxi  test  data  served  as  the 
basis  for  the  development  of  the  roll/pitch  control  system 
presented  in  this  paper.  Due  to  computational  problems  with 
the  air  cushion  model  of  the  computer  program,  the  controller 
designs  could  not  be  completely  verified;  but  expected  trends 
in  pitch,  roll  and  yaw  control  were  shown. 
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THE  DESIGN  OF  A TAKEOFF  STABILIZATION 
SYSTEM  FOR  AN  AIRCRAFT  WITH  AN  AIR  CUSHION  LANDING  SYSTEM 

Chapter  I 
Introduction 

In  the  low  speed  range  of  a takeoff  roll,  the  normal  air- 
craft controls  are  not  aerodynamically  effective;  hence,  the 
pilot  must  control  the  heading  by  differential  braking  and  wait 
until  the  ailerons  become  effective  to  control  roll.  During 
this  time,  the  landing  gear  dampens  most  of  the  pitch  and  roll 
oscillations  so  that  the  pilot  has  few  corrections  to  make  in 
the  latter  part  of  the  takeoff  roll.  However,  when  the  conven- 
tional landing  gear  Is  replaced  by  an  Air  Cushion  Takeoff  System 
(ACTS)  the  pitch  and  roll  damping  is  greatly  reduced.  This 
paper  will  use  the  Jlndivik  Remotely  Piloted  Vehicle  as  an  ex- 
ample of  an  air  cushion  vehicle  that  can  be  controlled  In  the 
low  speed  range  with  the  use  of  small  jet  thrusters  on  the  wing 
tips  and  a thrust  deflector  on  the  tail  section. 

The  Jlndivik  Remotely  Piloted  Vehicle  (RPV)  Is  an  Australian 
target  drone  that  can  be  launched  and  recovered  on  a runway.  At 
present,  the  takeoff  Is  accomplished  with  s takeoff  dolly,  as 
shown  In  Fig.  1,  that  provides  a wing  level,  attitude  and  direc- 
tional control.  At  lift  off  the  Jlndivik  separates  from  the 
dolly  and  the  dolly  brakes  to  a stop.  Recovery  of  the  Jlndivik 
Is  done  by  landing  on  a single,  four  Inch  wide  metal  skid  attach- 
ed to  the  fuselage.  Directional  control  during  landing  Is 
maintained  with  the  ailerons,  the  rolling  moment  thus  produced 
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Fig.  I.  Jlndlvik  on  a Takeoff  Dolly 


makes  Che  drone  ride  up  onto  an  edge  of  the  skid  and  turn  In 
the  direction  of  the  roll.  For  the  last  twelve  years  the 
Australian  Air  Force  has  used  the  Jindivik  in  this  configura- 
tion with  considerable  success. 

A Joint  project  by  Che  Australian  Air  Force  and  Che  United 
States  Air  Force  was  initiated  in  1972  to  incorporate  an  Air 
Cushion  Landing  System  on  the  Jindivik.  The  objectives  of  this 


project  were  to  convert  the  Jindivik  to  an  all-terrain  RFV  and 
to  advance  air  cushion  technology.  The  drone  end  air  cushion 
are  shown  in  Fig.  2.  Initial  low  speed  taxi  tests,  completed 
in  Australia,  show  chat  Che  aircraft  fitted  with  Che  air  cushion 
is  unstable  in  roll,  pitch,  and  directional  control  (yaw)  (Ref. 
13).  Therefore,  a Stability  Augmentation  System  (SAS)  will  have 
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Pig.  2.  Jlndivlk  and  Air  Cushion  Landing  System 
to  be  designed  and  incorporated  into  the  autopilot  before  the 
RPV  is  airworthy. 

Since  the  drone  was  designed  to  be  launched  from  a direc- 
tionally controlled  dolly.  It  was  not  designed  with  a rudder. 
Implementation  of  a rudder  during  this  project  would  require 
extensive  structural  changes  and  major  changes  to  the  autopilot 
and  ground  control  units.  Therefore,  a yaw  thruster  was  de- 
signed and  fitted  to  the  rear  of  the  fuselage  to  direct  the 
jet  exhaust,  thus  providing  a yawing  moment.  Roll  and  pitch 
control  will  be  provided  by  a vertical  roll  thruster  on  the 
front  tip  of  each  wing  pod.  The  roll  thrusters  will  be  ac- 
tivated alternately  to  control  rolling  moments  and  in  tandum 
to  control  pitch.  Since  the  roll  thrusters  are  on  the  tip 
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Pig.  3.  Block  Diagram  of  Autopilot  and  SAS  Control  Units 
of  the  wing  pods,  they  are  approximately  six  feet  ahead  of  the  centre 
of  gravity  of  the  drone  and  hence  can  produce  a moment  to  control 
the  pitch  attitude  to  some  degree.  As  well,  the  roll  thrusters 
can  be  directed  up  or  down  to  counteract  positive  or  negative 
rolling  and  pitching  moments. 

At  velocities  near  fifty  knots,  the  ailerons  and  flaps  be- 
come aerodynamically  effective  and  the  roll  thrusters  are  phased 
out  to  ensure  chat  the  vehicle  is  not  overcontrolled.  An  addi- 
tional advantage  of  turning  off  the  roll  thrusters  is  that  a 
smaller  gas  supply  is  required  for  the  thrusters.  Hence,  they 
can  be  used  with  a gas  bottle  rather  than  bleed  air  from  the 
engine.  This  arrangement  will  greatly  reduce  the  airframe  and 
engine  modifications  required  for  implementation. 
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i The  stability  augmentation  system  will  be  designed  to  use 

I 

f 

I the  existing  sensors  In  the  autopilot  to  control  the  roll  and 

pitch  thrusters.  The  SAS  unit  will  be  placed  In  the  feedback 
control  loop  between  the  autopilot  and  actuators,  as  shown  In 
Fig.  3. 

i 

This  thesis  Is  organized  In  the  following  manner:  Chapter 
II  develops  the  equations  of  motion  and  aerodynamic  stability 
derivatives.  Chapter  III  describes  the  air  cushion  model,  Chap- 
ter IV  discusses  the  controller  design.  Chapter  V contains  a 
description  of  the  computer  program  and  the  simulation  results, 
and  Chapter  VI  the  conclusions  and  recommendations. 


The  following  six  simultaneous  non-llnear  differential 
equations  fully  describe  the  motion  of  the  Jlndlvlk  RPV.  The 
positive  sense  of  the  variables  Is  in  the  direction  of  the  arrows 
in  Fig.  4. 

~ VB. +Fax 

+UR-WP)  ^ +Fy£xr 

^('vV  -UQ+VP:)  = cos^i 

X/K/x  P-Xx^R'^x^PQ.  -♦•(X^«*”XypRQ  = ■*‘^£xr  (^) 

Xyy<!^-P^xx”’Izs)PR'Flx^(P*'-"R^)  <5) 


R "Xx^F  •t-OCy^-rxx)PQ.'''^X^QR=  '^‘^eXT  (6) 

The  equations  are  first  order  in  U,  V,  W.,  P,  Q,  and  R with  the 
added  kinematic  relationships. 

p = 4>  - V © (7) 

Q =:  ^ COS  + y COS  o sr<A/ 4*  (8) 

= 9<^os  e cos  4>  - e sm  4>  O) 


The  assumptions  used  In  the  derivation  of  these  equations  were: 
(1)  the  aircraft  Is  a rigid  body,  (2)  the  mass  of  the  aircraft 
Is  constant  for  the  duration  of  the  analysis,  (3)  gravity  Is 
constant,  (4)  the  earth  Is  an  Inertial  reference,  (5)  and  there 
is  body  axis  symmetry  about  the  x-z  plane  (i.e.  ,Ix/*0=Iy2  )• 

Equations  (1)  to  (9)  are  Included  in  the  subroutines  of  the 
"EASY  Dynamic  Analysis  Computer  Program  to  Aircraft  Modeling" 
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Fig.  4.  Definitions  of  Vector  Components  In  the 
Equations  of  Motion 

(Ref.  5)  which  was  used  to  simulate  the  takeoff  motions  of  the 
Jlndlvlk.  The  Inputs  required  by  EASY  are  the  mass,  inertias, 
geometry  of  the  aircraft,  and  tne  aerodynamic  and  external  forces 
acting  on  the  drone.  The  air  cushion  system  Is  considered  to  be 
a prime  generator  of  the  external  forces  ana  moments  (aside  from 
engine  forces  and  moments)  and  Is  described  in  the  next  chapter. 
The  mass,  inertias  and  geometry  are  readily  available  from  air- 
craft blue  prints  and  reference  manuals;  and  the  aerodynamic 
forces  and  moments  can  be  computed  from  wind  tunnel  model  data 
and  theoretical  methods. 

The  aerodynamic  forces  and  moments  can  be  written  as: 

L = ^ S 

p =c.t,  fS 
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where 


Fy-  CyfS 

(12) 

(13) 

??= 

(14) 

7fl-  C/*i  ‘j'S 

(15) 

Cl- 

(16) 

Ct)= 

(17) 

Cy= 

(18) 

(19) 

Cm.  - sv,) 

(20) 

Cyjrv  * 

(21) 

nts  Cl  , Ct>  , Cy  , C4  , C.^  , and 

^/»rv 

are  non-dioenslonal.  By  determining  the  coefficients  for  every 
flight  condition,  the  aerodynamic  forces  and  moments  can  be  cal- 
culated and  added  to  the  external  forces  and  moments  to  produce 
the  aircraft  motions.  These  calculations  are  done  by  the  EASY 
program  but  the  program  requires  all  the  aerodynamic  stability 
derivatives  (all  the  functional  relationships  which  determine 
the  force  and  moment  coefficients).  The  remainder  of  this  chap- 
ter will  deal  with  the  derivation  of  the  stability  derivatives. 
These  derivatives  arc  derived  in  the  stability  axis  system  as 
defined  by  Blakelock  (Ref.  2). 
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Reference  13  contains  wind  tunnel  data  for  the  Jlndlvlk  In 
various  configurations,  Including  one  when  fitted  with  the  Air 
Cushion  Recovery  System  (ACRS).  The  ACRS  Is  the  air  cushion 
trunk  with  which  the  drone  lands,  but  it  also  has  an  Air  Cushion 
Takeoff  System  (ACTS)  trunk  with  which  It  takes  off.  After 
takeoff,  the  ACTS  Is  disengaged  and  drops  to  the  ground.  Both 
trunks  are  the  same  shape  with  the  ACTS  being  about  21Z  larger 
In  all  dimensions.  Since  no  wind  tunnel  data  was  available  for 
the  ACTS,  the  data  for  the  ACRS  was  extrapolated  by  percentages 
and  assumed  to  be  fairly  accurate  for  the  ACTS.  An  example  of 
the  estimation  technique  Is  that  the  increase  in  the  frontal 
area  of  the  aircraft  due  to  the  replacement  of  the  ACRS  by  the 
ACTS  was  6Z:  therefore,  the  values  of  the  coefficient  of  drag 
( C.U  ) were  Increased  by  7Z.  Since  the  trunk  does  not  generate 
lift,  the  coefficient  of  lift  ( ) was  not  affected,  nor  was 

the  coefficient  of  side  force  ( Cy  );  symmetry  about  the  x-z 
axis  meant  that  the  coefficient  of  yawing  moment  ( ) was 

not  affected.  The  pitching  moment  coefficient  (C^  ) and 
rolling  moment  coefficient  ( ) were  affected  by  the  per- 

centage that  the  increased  drag  affected  those  moments.  Thus, 

Co  . , Cy  . Cji  , , and  C/fy  can  be  empirically  de- 

termined as  functions  of  the  angle  of  attack  ( ><  ) t the  sideslip 
( ^ ),  and  the  elevator  deflection  ( fe  ) • I"  other  words  , 

^ ■ 4?f  • ^ ^ ^ 

, and  can  be  found  from  the  wind  tunnel  data. 

Non-dimensional  derivatives  were  calculated  because  they 
provided  a means  of  checking  typical  values  and  signs  with 
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Roekam  (Ref.  17)  and  Blakelock  (Ref.  2).  Before  entering  Che 
derivatives  into  the  EASY  program,  they  were  dimenslonallzed. 

At  low  airspeeds  the  heave  motion  of  the  air  cushion  can 
create  angles  of  attack  beyond  the  stall  limit,  but  at  these 
speeds  aerodynamic  contributions  to  the  aircraft  dynamics  are 
small. 

Stability  Derivative  Derivation 

Cl. 

From  a curve  of  C^vs.  ^ of  the  wind  tunnel  data  it  can 
be  shown  that 

Cl  = C^^+.^  ««  (2 

andC^^.^^  can  be  determined  directly  as  the  inter- 
cept and  slope  of  the  curve. 


.Q.EL 


From  a drag  polar  of  vs.  C-q  ft  can  be  shown  chat 

Co  - KCj-  (23) 

where  and  K are  determined  by  a curve  fit  of  wind  tunnel 


data  of  C-0  and 

C(,  . Substituting  for  Ci_ 

~ Cq^ -h -h  Ci.,^ 

= Cd«+ K (£,>  ac..c  2 -. + <') 

(24) 

differentiating 

•^=  K «■) 

(25) 

= aKCL<Cc..<,+CL,-<) 

so 

(26) 

Roskam  (Ref.  17 

, pgs  4.12,  1.18,  4.25)  shows  that  for  velo.clcles 

below  300  ft/sec  the  variation  of  lift,  pitch  moment,  and  drag 
with  velocity  is  zero.  Thus, 


(27) 


APIT/GE/EE/77D-43 
The  quanCltles 

SCo  , ^ SCi.  , <^C/m,  , , <^C/>n 

d S'c  d S9  6%q,  d S9  d •<  dS^e  "dTf" 

were  all  computer  from  curve-fits  of  the  wind  tunnel  data. 

From  Roskam  (Ref.  17)  it  can  be  8ho\m  that  the  angle  of 
attack  of  the  tail  in  dovnwash  is 

— oi  — L ^ 

and  for  a particular  angle  of  attack 


(28) 


a 


(29) 


now  the  change  in  lift  coefficient  on  the  tail  due  to  dovnwash  is 


The  change  in  aircraft  lift  is 

s 


c.  = 

/)R  C.OS  A Cl-c^-d 

where  ^ *4  tr 

^ Tr/»R  — TT 

= 5.73 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


X is  sweep  angle  and  /)/^  is  aspecf  ratio  (Ref,  7) 
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The  wing  contribution  to  Is  considerable  but  can  not  be 

estimated  (using  Roskam,  DATCOM,  etc.)*  So  a "typical"  value 
from  Roskam  of  /•S'rnA  was  used;  fortunately  this  deriva- 

tive Is  of  minor  Importance  (Ref.  17,  p.  4 >114). 

Cdn  s 

The  contribution  of  the  wing  was  neglected  because  It  will 
be  negligible  with  respect  to  the  tall  contribution.  The  correc- 
tion to  the  pitching  moment  due  to  downwash  on  the  tall  Is 


and 


and 


a<  c u S 

S-t 

C U S 


(36) 

(37) 

(38) 


^ changes  the  angle  of  attack  on  the  tall  by,g^  radians 
(for  quaslstatlc  conditions) 


so 


and 


M 

u 


(39) 

(40) 

(41) 


differentiating  Ssii 

Q^>a.u  ^ s u 

and  the  contribution  of  the  wing  body  Is  negligible  In  comparison 
to  the  tail  (Ref.  17,  p.  153) 


now 


Cta-  4 aji  ^ = a.a^t  s.- 
V c as 


(42) 
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I 

j 

^ ) 

r t 

f I 


C'MfX 


The  Domenc  on  the  tall  Is 


7n^ 


(A3) 


and  the  change  In  moment  due  to  a change  In  angle  of  attack  Is 

acs^ 

v<^f>VAic  acs 

since  the  wing  contribution  Is  negligible  with  respect  to  the 
tall  (Ref.  17,  p.  153) 


(AA) 

(A5) 

(A6) 


i 


4 ^ 


C^S 


(A7) 


Derivatives 


Wind  tunnel  data  gaveOy  , , and  ^/n  vs.  ^ for  values 

of  / p/i7*,  but  in  any  takeoff  with  crosswind  the  sideslip  will 
normally  exceed  7°.  The  normal  takeoff  procedure  will  be  to 
initially  line  the  aircraft  Into  the  relative  wind  at  the  center- 
line  and  change  the  heading  as  the  aircraft  gains  speed.  This 
procedure  should  keep  ^ within  +30®  and  the  present  data  can  be 
curve  fitted  and  extrapolated  to  this  value.  Consequently,  ex- 
pressions can  be  obtained  for  t i <ui<)  from  the  data. 

The  ^ derivatives  have  been  assumed  to  be  .zero  (Ref.  17).  For 
the  ^ , p , and  A.  derivatives  the  effect  of  sidewash  on  the  tall 
has  been  neglected. 
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where  ^ ^ is  Che  mean  height  of  Che  fin,  and  the  effect  of  sldewash 
has  been  neglected 


so 


= TK^Ji^ 


(5o: 


The  wlag  contribution  Is  in  two  parts  due  to  lift  and  drag.  For  positive 
p>  the  angle  of  attack  is  Increased  on  the  right  wing  and  decreased  on 
Che  left;  thus  Inclining  and  changing  Che  lift  and  drag  vectors  of  each 


wing  section.  The  inclination  angle  is  ^p=-P  Y where  w is  the  spanwlse 

u ^ 


coordinate  of  the  section.  The  change  in  lift  is 

Jls 


(51 


and 


pv  Cci' 


(52 


ire  C</y 


vliere  Oay  is  the  area  of  the  wing  section- 
So  the  change  in  the  yawing  moment  due  to  lift  is 


= -y 


(53 


i 


and  the  total  yawing  moment  of  the  wing  Is 


p*"  c::. 


= — ^ ^ 

3^  u^ 


=.  - 

dp  I ^ U.^ 

^ '\  y /,  /)  ^ -2J{.  - _ ^f-ac  p fc>^c 


the  change  In  drag  Is 


AOjupr  ^ -% 

- - s ^ 


the  corresponding  change  In  yawing  mcment  Is 


A - 

^ T^Aut^T^  ~ 


5 

^hr- 


p_y<^y 


S ££iU2-y-‘="'y 
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Cjip. 

C-yple  often  negligible  (Ref.  17,  p.  170)  and  the  tall  la  the  inajor 
contributor.  Let  the  mean  change  In  angle  of  attack  of  the  vertical 
stabilizer  (fin)  be 


« —-EAE: 
u 

where  la  the  mean  height  of  the  fin 
now  ^ Cy  ^ ^ fz  ^ 


so  the  change  on  the  side  force  coefficient  of  the  aircraft  ls.(sldewash 
la  neglected) 

SO  X - Of  ^<~^F  ( 

dp  s u 


Cyo=  ^ r ~ 

b Op  < L 
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The  wing  is  the  only  major  contributor  (Ref.  17,  p.  170).  Etkln 
(Ref.  9)  shows  that- the  rolling  velocity,  p , produces  a change  In  the 
angle  of  attack  of  each  wing  section  which  is  proportional  to  the  span, 
l.e., 

U (73) 

where  y la  the  spanwlse  coordinate  of  the  wing  section.  Then  the  lift 
distribution  on  the  wing  due  to  rolling  Is  estimated  to  be 

y 

Etkln  (Ref.  9)  changes  the  triangular  lift  alstrlbution  to  a 
sinusoidal  distribution  to  account  for  the  loss  of  lift  at  the 
wing  tips  due  to  spanwlse  flow  around  the  wing  tips.  However,  since 
the  Jlndlvlk  has  large  tip  tanks  the  lateral  airflow  will  be  mlnimltiid 
and  the  lift  distribution  will  be  closer  to  the  triangular  distribution 
shown. 

The  change  in  lift  will  be 

AL  ^ 2.  f 

and  the  change  In  rolling  moment  due  to 

17 
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y 


<^py 

u 

this  lift  will  be 

c py  c^y 
U 


(74) 


(75) 


Will 

t 


( 

1 
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thus 

(76) 

so 

d p 2.^  u 

(77) 

and 

<X!c 

(78) 

so 

^ Cc^ 

^ /as  la 

(79) 

Derivatives 


The  call  is  Che  prime  contributor  (Ref.  17,  p.  174).  The  change  In 
fin  angle  of  attack  due  to  a yaw  rate  is 


^o<.  - -A, 


(80) 


so  the  change  In  aide  force  due  to  the  Call  Is 

7 “ 


u 


4f^ 

c)/c  ^ S (( 


(81) 

(82) 

(83) 


Contributions  are  from  the  wing  and  tall.  The  side  force  on  the  .all 
acts  at^/r,  Che  mean  height  of  Che  fin 
and  ^ ^ ^ 


^ s a 


(84) 
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now 


differentiating 


= Sr: 

S (X 


^ ‘ S W 


(85) 


(86) 


and 


— L_  c^Sf 

b Sfb 

^ ^ ;^>p 

Sb^ 

A positive  A.  also  increases  lift  on  the  left  wing  and  decreases  it  on 
the  right;  the  change  in  lift  on  each  wing  section  is 
^ -oa,etV>o  = -£>  ^ 

_ i ^ (^uf  Cl  c 

^ -X  ^ (^y)  ^ 

now  the  change  in  rolling  moment  due  to  two  sections  at  y is 


(87) 


(88) 


(89) 


(90) 


and  the  total  rolling  moment  change  for  the  wing  is 

- p/t^  Cl  c b'*' 


- ^ c 


3ZU^ 


(91) 
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. a^Cct"' 


CA  4j iX 


1 6 a' 
I C>u'^ 


(tjiA  . = M. ^ ~ ■^Cf.h 
^ b at  a a 


summing  the  components 


dc  = O.^LcLf=-  Ct.i 

BU 


^ILO. 


The  tail  and  wing  contribute  toO»/j.  Knowing  that  the  change  in  the 


fin  angle  of  attack  is 


^ p. 

u 

and  the  noment  arm  of  the  tall  is  then 


■ g> 

5 ^ 
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A positive  /t  Increases  the  drag  on  the  left  wing  and  decreases  drag 
on  the  right  wing.  The  change  In  drag  on  eac\  wing  section  Is 

-X.-p  Cj,  cAy 

=:X  y^cc/y 

and 


SO  the  change  In  yawing  moment  Is  i 

^ 3 

A9f  y 


ClOO) 

(101) 


P C-o  ^ 

(102) 

now 

- — p - - 

^ Cjj  /I  c l>^ 

(103) 

3-^  3^ 

ife  uy 

and 

^ i)  5^b 

C-[>yx.  ii 

(104) 

■> 

8 U 

summing  the 

components 

• 

C/n/»  ~ 

(105) 

su 

Once  all  the  derivatives  had  been  determined  or  estimated.  Ref.  15  was 
used  to  convert  to  dlmentlonal  body-axis  derivatives.  The  derivatives 
were  the  written  Into  the  EASY  program. 
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Chapter  111 
The  Air  Cushion  Model 

The  air  cushion  model  used  in  this  analysis  Is  a truncated  version 
of  an  ACLS  model  that  was  designed  by  Foster-Miller  Associates  Inc.  of 
Waltham,  Massachusetts  for  the  National  Aeronautics  and  Space  Adminis- 
tration (NASA)  (Ref.  4). 

The  basic  ACTS  configuration  Is  shown  In  Fig.  5.  The  model  in- 
cludes four  primary  subsystems:  (1)  the  fan,  (2)  the  feeding  system, 

(3)  the  trunk,  and  (5)  the  cushion.  Air  from  the  fan  flows  through 
the  ducts  and  plenum  (feeding  system)  and  eaters  the  trunk.  The  trunk 
has  several  rows  of  orifices  that  exhaust  both  to  the  cushion  and  the 
atmosphere.  Thus,  the  airflow  from  the  trunk  has  two  components, 
one  entering  the  cushion  and  the  other  leading  directly  to  the  atmos- 
phere. The  cushion  flow  exhausts  to  the  atmosphere  through  the  clearance 
gap  formed  between  the  trunk  and  ground.  In  addition  to  the  basic  flows 
described  above,  two  other  flows  have  been  included  in  the  model.  These 
are  the  plenum  bleed  flow  and  the  direct  cushion  flow.  Plenum  bleeding 
causes  some  of  the  air  to  flow  directly  from  plenum  to  atmosphere,  and 
has  been  used  In  some  designs  to  improve  the  dynamic  characteristics  of 
the  air  supply  system.  Direct  flow  from  the  plenum  to  the  cushion  can 
also  Improve  dynamic  response.  A pressure  relief  valve  Is  also  Included 
In  the  basic  configuration.  It  allows  additional  flow  to  vent  from  the 
plenum  whenever  the  pressure  exceeds  a preset  level , and  thus  Improves 
stability  by  reducing  fan  stall. 

The  support  force  acting  on  the  aircraft  Is  made  up  of  two  components. 
The  first  occurs  due  to  the  cushion  pressure  acting  over  the  cushion  area. 
The  second,  which  cornea  about  only  during  ground  contact,  is  given  by  the 
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Flow 

From 

Fan* 

Bleed  Flow 


Fig.  5.  Basic  ACTS  Configuration 
contact  pressure  acting  over  the  trunk  contact  area.  The  support 
force,  In'  general,  also  gives  rise  to  a moment,  given  by  the  pro- 
duct of  the  force  and  Its  distance  from  the  CG  of  the  aircraft. 

In  plan,  the  cushion  has  an  oval  shape,  made  up  of  a rec- 
tangular section  with  semicircular  ends.  The  lengths  a and  b are 
the  horizontal  and  vertical  spacing  between  the  points  of  attachment 
of  the  trunk  to  the  aircraft  body.  The  initial  (undeformed)  trunk 
shape  Is  defined  by  the  above  two  parameters  and  the  perimeter 
and  height  Hy  as  shown.  is  the  (uniform)  spacing  between  the 
rows  of  peripherally  distributed  orifices.  The  number  of  the 
orifices  Is  selected  independently  by  the  number  of  orifice  rows 
Nf  and  the  number  of  orifices  per  row  The  cushion  volume  consists 

of  two  parts:  an  active  (dynamically  varying)  region  and  a dead  (static) 
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region.  The  active  volume  depends  on  the  crunk  shape  and  ground  pro- 
file. The  dead  volume,  which  Is  a design  variable,  Includes  recesses 
In  the  cushion  cavity  as  shown. 

The  forces  transferred  to  the  aircraft  act  through  the  cushion 
and  crunk.  To  help  calculate  these  forces,  the  crunk  and  cushion  are 
divided  Into  segments  as  shown  In  Fig.  6.  Each  straight  section  of  the 
cushion  and  trunk  Is  divided  Into  M rectangular  segments,  while  each 
curved  end  Is  divided  into  N pie-shaped  segments.  Thus,  Che  total 
number  of  segments  Is  2 (M  + N) . All  cushion  and  trunk  parameters 
are  calculated  first  for  each  segment  and  then  summed  to  give  their 
total  system  values. 

The  dynamic  analysis  of  the  vehicle  system  is  best  derived  with 
the  help  of  two  orthogonal  coordinate  frames  of  reference:  a coordi- 
nate frame  fixed  In  space  (inertial  frame) , and  a coordinate  frame 
fixed  to  the  vehicle  (vehicle  frame)  with  origin  at  Che  aircraft  CG. 

The  reason  for  two  frame's  can  be  appreciated  by  recognizing  thst; 

(a)  Newton's  law  for  translation  motion  requires  that  the 
CG  acceleration  be  expressed  relative  to  the  inertial  - 
frame. 

(b)  The  corresponding  law  for  rotational  motion,  while  valid 
In  both  Inertial  and  vehicle  frames.  Is  applied  more  con- 
veniently In  Che  vehicle  frame,  because  rotational  Inertia 
about  any  vehicle  axis  Is  constant,  while  the  rotational 
Inertia  about  any  Inertial  (fixed)  axis  varies  with  air- 
craft position. 

Accordingly,  the  two  frames  of  reference  have  been  defined  as 
shown  in  Fig.  7.  The  vehicle  frame  with  origin  at  the  aircraft  CG 
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Fig.  6.  Division  of  Trunk  Into  Segments 
has  roll,  yaw  and  pitch. axes  x,  y and  z,  respectively,  fixed  to  the 
Aircraft  body  as  shown.  The  inertial  frame  has  corresponding  axes 
X,  Y and  Z fixed  in  space.  The  two  frames  coincide  only  when  th^  air- 
craft has  not  undergone  any  rotation  from  equilibrium. 

In  the  analysis,  the  actual  runway  profile  underneath  the  ACTS 
is  approximated  by  segments  that  coincide  in  plan  with  those  of  the 
trunk  and  are  parallel  to  the  cushion  hard  surface  as  shown  in  Fig.  7. 
With  this  model,  all  pressure  forces  act  parallel  to  the  vehicle  yaw 
axis  so  that  the  segment  torque  components  about  the  aircraft  CG  can 
be  easily  computed  by  multiplying  the  segment  force  by  the  fore-and- 
aft  and/or  lateral  separation  between  the  segment  and  the  CG. 
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Fig,  7.  Inertial  and  Vehicle  Coordinate  Frames 
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The  analytical  model  of  the  ACTS  consists  of  a set  of  equations 
which  when  solved  determines  the  pressures,  flows,  forces,  and  motion 
of  the  system  for  various  aircraft  and  runway  parameters.  The  overall 
ACTS  system  is  divided  Into  two  interrelated  systems:  the  flow  system 
and  the  force  system.  These  systems  are  shovm  in  Fig.  8 and  Fig.  9. 

The  flow  system  establishes  the  pressure-flow  relationship  for  various 
subsystems  of  the  ACTS.  The  force  system  establishes  the  corresponding 
force-motion  relationships.  The  interdependence  of  the  two  systems  comes 
about  because  the  crunk  deflection  obtained  from  the  force  system  changes 
the  volumes  and  orifice  areas  that  form  part  of  Che  flow  system.  Simi- 
larly, the  cushion  and  trunk  pressures  found  from  the  flow  system  give 
rise  to  forces  and  moments  that  form  Inputs  to  the  force  system. 

The  Trunk  Model 

The  major  component  of  the  ACTS  model  is  the  trunk  model  because 
It  determines  the  trunk  shape  parameters  (volume,  and  orifice  and  con- 
tact areas) , contact  pressure  distribution  and  damping  that  form  Inputs 
to  the  ACLS  flow  and  force  systems. 

Trunk  Shape.  In  past  work,  two  analytical  models  have  been  developed 
for  the  trunk  shape:  the  Membrane  Trunk  Model  (Ref.  8)  and  the  Frozen 
Trunk  Model  (Ref.  6).  The  shortcoming  of  both  these  analyses  was  that 
they  modeled  the  side  and  end  segments  of  the  trunk  In  the  same  way 
while  test  data  now  confirm  Chat  the  shorter  curved  end  segments  (front 
and  rear)  behave  very  differently  from  the  longer,  straight  side  seg- 
■a  nta.  F.ig.  10  shows  the  trunk  cross  section  measured  at  the  center 
of  the  side  and  end  segments  as  the  load  on  Che  cushion  Is  Increased. 

The  entire  side  segment  tends  to  bow  outward  and  avoid  ground  contact, 

I 

while  the  end  segment  remains  virtually  fixed,  except  for  a flattening 
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Fig.  9.  ACTS  Force  System 

in  the  region  that  actually  touches  the  ground.  This  difference  in 
behavior  occurs  because  the  front  segment  is  much  smaller  than  the 
side  segment  and  Is  curved.  When  the  cushion  pressure  increases  due 
to  an  Increase  in  the  load,  the  radially  outward  force  causes  the  oval 
trunk  planform  to  become  more  circular,  as  shown  in  Fig.  11.  This 
causes  a hoop  tension  force,  T,  to  act  around  the  trunk  periphery. 

In  the  side  segments,  this  force  acts  substantially  normal  to  the 
side  excursion,  ^g,  so  that  Its  component  resisting  the  motion  is 
negligible  and  the  side  segment  can  chus  bow  outwards  relatively  un- 
restrained. In  the  end  segments  the  situation  Is  different,  since 
the  curvature  of  the  segment  causes  the  hoop  tension  to  have  a much 
higher  component  opposing  the  motion  so  that  outward  motion  of  the 
trunk  ends  Is  much  smaller. 
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Fig.  11.  Outward  Excursion  of  Trunk  Segments 
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Since  hoop  tension  has  very  little  effect  on  side  trunk  notion » 
the  side  segnents  can  be  considered  as  simple,  two-dimensional  mem- 
branes, as  done  In  the  Membrane  Trunk  Model.  On  the  ocher  hand,  the 
fact  that  hoop  tension  restrains  ("freezes")  the  trunk  ends  suggests 
that  these  segments  be  modeled  by  the  Frozen  Trunk  Model.  Thus,  the 
logical  step  in  trunk  model  Improvement  Is  to  combine  the  two  exist- 
ing models  and  form  Che  Hybrid  Trunk  Model,  in  which  the  sides  are 
represented  by  the  Membrane  Model  and  the  ends  by  the  Frozen  Model. 

The  Hybrid  Trunk  Model  is  essentially  a limiting  case  analysis 
6£  trunk  deflection.  In  general,  best  results  will  be  obtained  at 
the  middle  of  the  respective  segments,  i.e.,  at  the  center  of  the 
side  segments,  where  the  crunk  behaves  very  much  like  an  ideal  mem- 
brane, and  at  the  center  of  the  end  segments,  where  the  trunk  shape 
Is  truly  fixed.  In  the  transition  region  (at  and  near  where  the 
segments  meet) , Che  trunk  will  exhibit  properties  of  both  the  mem- 
brane and  frozen  trunk  approximations. 

Contact  Pressure.  In  addition  to  trunk  and  cushion  shape,  the 
trunk  model  also  determines  the  pressure  distribution  in  the  ground 
contact  zone.  The  analysis  for  pressure  distribution  Is  complicated 
by  Che  fact  that  two  separate  effects  must  be  considered:  direct 
trunk-ground  contact  caused  by  the  trunk  pressure  forcing  the  trunk 
against  the  ground,  and  airflow  through  the  trunk  holes  Into  Che 
interstices  that  remain  In  Che  contact  zone. 

When  two  bodies  In  contact  are  acted  upon  by  a force,  F,  the 
actual  contact  occurs  at  a number  of  discrete  regions  rather  than 
over  the  whole  area,  due  to  the  Inherent  roughness  of  the  contact- 
ing surfaces  as  shown  In  Fig.  12.  Because  the  number  of  contact 
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regions  Is  large  for  a "snooth"  surface,  it  was  convenient  to  define  an 
average  contact  pressure,  • F/A,  acting  as  though  the  bodies  were 
touching  uniformly  over  the  entire  area,  A.  In  fact,  P^y  equals  the 
trunk  pressure,  P^.  For  purposes  of  trunk  outflow  calculation,  the 
pressure  profile  in  the  non-contacting  regions  is  approximated  by  a 
linearly  decreasing  relationship  as  shown  in  Fig.  12.  The  driving 
pressure  for  flow  through  any  trunk  hole  is  rhus  g .'>en  by  the  differ- 
ence between  the  trunk  pressure  and  the  gap  rr  .f.  • 't  that  location. 
This  pressure  distribution  model  has  been  compared  to  experimental 
dati;  and  has  t4en  quite  accurate  (within  10%)  (Ref.  4). 

Trunk  Damping.  In  dynamic  operation,  the  trunk  is  deformed  cyclically 
both  in  tension  and  flexure,  and  energy  dissipation  in  the  trunk  material 
gives  rise  to  a damping  force  which  opposes  the  strain  :<-  te.  Because 
the  present  trunk  analysis  does  not  solve  for  strain  (auJ  hence  strain 
rate) , a damping  model  that  links  trunk  material  properties  directly 
to  trunk  damping  forces  cannot  be  developed.  An  alternate  approach 
in  which  the  damping  characteristics  are  modeled  by  dimensional  analysis 
(similarity)  based  on  test  data  thus  appears  more  appropriate.  In 
keeping  with  the  method  of  approach  outlined  earlier,  the  trunk  is 
divided  into  segments  (Fig.  13)  and  a series  of  dashpots,  one  for  each 
segment,  is  included  in  the  model  such  that  the  segment  damping  force 
is  proportional  to  the  vertical  velocity  of  the  trunk  segment. 

Each  dashpot  models  the  energy  dissipation  characteristic  of  the 
trunk  segment.  Although  all  parts  of  the  trunk  dissipate  energy,  the 
major  contributions  will  come  from  those  parts  that  undergo  high 
stress  reversals,  since  the  strain  rate  is  .highest  in  these  sections. 
Observations  of  a trunk  in  dynamic  operation  suggest  that  the  high 
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stress  reversal  regions  lie  along  the  periphery  of  the  trunk-ground 
contact  zone,  because  it  is  here  Chat  the  rate  of  change  of  crunk 
slope  (and  hence  stress)  Is  high  and  constantly  changing  with  the 
time  as  Che  contact  area  changes.  As  a first  order  approximation, 
the  damping  model  derived  here  assumes  Chat  all  the  energy  dissipa- 
tion In  Che  trunk  Is  concentrated  along  the  Crunk-ground  contact 
periphery  so  that  Che  damping  coefficient  of  each  dashpot  depends 
on  the  perimeter  of  the  ground  contact  zone.  This  means  that  when 
a segment  is  not  contacting  Che  ground  it  has  zero  damping  and  when 
it  is  contacting  the  ground  it  has  a damping  coefficient  proportional 
to  the  contact  perimeter. 

Model  Synopsis 
The  Flow  System 

(a)  The  fan  is  characterized  by  a static  pressure  rise  element  for 
forward  and  back  flow  in  series  with  an  Inertance  (duct)  and  a 
capacitance  (volume). 

(b)  The  trunk  and  cushion  volume  are  found  from  the  Hybrid  Trunk  Model, 
which  characterizes  Che  side  trunk  segment  as  an  ideal  two-dimen- 
sional membrane  and  the  end  segment  as  a "frozen"  trunk. 

(c)  The  orifice  areas  between  the  trunk  and  cushion,  trunk  and  atmos- 
phere and  cushion  and  atmosphere  are  found  from  the  crunk  shape 
as  predicted  by  the  Hybrid  Trunk  Model,,  along  with  the  cushion 
orientation  and  ground  profile. 

(d)  The  pressure  within  Che  cushion,  trunk  and  plenum  is  considered  to 
be  uniform. 

(e)  The  pressure  In  the  trunk/ground  contact  zone  is  found  from  the 
triangular  profile  given  by  the  Hybrid  Trunk  Model. 

36 


t 


US' 


AFIT/GE/EE/77D-43 

(f)  The  flow  through  the  plenum,  trunk  and  cushion  is  governed  by 
Che  unsteady  state  flow  continuity  equation  In  which  the  air  is 
assumed  to  behave  like  a perfect  gas  and  follow  a polytropic  ex- 
pansion relationship. 

(g)  The  flow  through  all  orifices  is  found  from  the  incompressible 
flow  square-law  orifice  equation. 

Die  Force  System 

(a)  The  mean  contact  pressure  in  the  trunk/ground  contact  zone  is 
equal  to  the  trunk  pressure. 

(b)  The  crunk  contact  area  and  location  relative  to  the  aircraft  CG 
is  found  from  the  trunk  shape  predicted  by  the  Hybrid  Trunk  Model. 

(c)  The  cushion  area  and  location  relative  to  the  aircraft  CG  is  found 
from  the  Hybrid  Trunk  Model.  In  width,  the  cushion  extends  between 
the  lowest  (ground  tangent)  points  of  the  side  crunk  segments.  In 
length,  it  extends  between  the  ground  tangent  points  of  Che  end  trunk 
segments,  or,  if  in  ground  contact,  between  the  inner  edges  of  the 
contact  zone. 

(d)  The  total  forces  and  moments  acting  on  the  aircraft  occur  due  to 
Che  mean  Crunk  contact  pressure  acting  over  Che  contact  area,  Che 
cushion  pressure  acting  over  the  cushion  area,  aerodynamic  drag 
and  trunk  damping  losses  caused  by  aircraft  heave  motion,  and 
trunk-ground  friction. 

(c)  The  forces  and  moments  are  found  by  dividing  the  cushion  (and  crunk) 
into  segments,  approximating  the  actual  ground  profile  underneath 
the  cushion  by  a similar  set  of  segments  parallel  to  the  cushion, 
computing  the  cushion  and  contact  pressure  forces  and  moments  for 
each  segment,  and  then  summing  them  Co  dec«rmine  Che  total  force 


i 


In  the  vertical  direction  to  the  aircraft  CG. 

(g)  Angular  accelerations  in  pitch  and  roll  are  obtained  by  applying 
the  theorem  of  moment  of  momentum  about  the  aircraft  pitch  and  roll 

axes. 

(h)  A coordinate  transformation  la  carried  out  to  express  vehicle  frame 
velocities  and  accelerations  in  terms  of  Euler  angles  and  their  de> 
rlvatives. 

(i)  Tlie  moment  of  momentum  equations,  expressed  in  terms  of  Euler  angles 
arc  Integrated  to  give  the  angular  position  of  the  aircraft  as  a 
function  of  time. 
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Chapter  IV 
Controller  Design 

During  the  low  speed  portion  of  the  takeoff  roll  (to  approxi- 
mately 50  knots) , the  controls  available  are  a yaw  thruster  on  the 
rear  fuselage  and  vertical  roll  thrusters  on  each  wing  dp.  The 
roll  thrusters  can  be  directed  up  or  down.  During  the  takeoff 
sequence,  the  most  unstable  mode  of  the  aircraft  Is  the  roll  mode. 

Therefore,  it  was  decided  to  control  this  mode  and  observe  the 

control  that  was  applied  to  the  pitch  mode  through  the  inertial  * 

cross-coupling.  Also,  during  takeoff  the  yaw  angle  will  be  con- 
trolled by  the  yaw  thruster  on  the  rear  of  the  fuselage.  With 
roll  and  yaw  controlled,  Eqn  (4),  which  Is  rewritten  here,  can 
be  simplified. 

I /x/x  P -hPa)  + 

Controlling  roll  and  yaw  means  that  R will  be  small  and  and 
8(^  (products  of  small  numbers)  will  be  small.  Cl  can  be  consid- 
ered small  because  the  takeoff  starts  with  zero  Initial  conditions  ' 

on  P , , and  P . With  the  above  simplification  and  assuming 

that  any  roll  Inputs  from  the  ground  profile  are  Impulsive,  then 
the  roll  moment  generated  by  the  roll  thruste''s  can  be  written  as 
Xxx  ^ ^ 


II 

(107) 

so 

=CFTC-t') 

(108) 

where 

C = 4^  = 8-18  X IO-’ 

•**xx 
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now  let 

/)(,=  4* 

(109) 

• 

/<x.-  ^ 

(110) 

and 

LL(t)- 

(111) 

so  Eqn  (108) 

can  be  rewritten  as 

- 

Jijt) 

(112) 

and 


so 


from  Eqn  (110)  (113) 


5 3 '■5] 


(114) 


^ e ^ 

For  a minimum  time  Performance  Index  and  for 


an 


optimal  control  for  this  system  can  be  shovm  to  be  a bang>bang 
control  (Ref.  14,  pgs  245-248).  In  other  words,  the  control  is 
a maximum  (either  positive  or  negative)  whenever  It  Is  applied. 
Since  the  eigenvalues  for  ^ are  both  zero.  Theorems  5" 4-1  and 
5*4-3  of  Kirk  (Ref.  14)  show  that  an  optimal  control  exists.  Is 
unique  and  has  at  most  cne  switching. 

Therefore,  the  control  for  a specified  initial  state  must  be 

LL  0^—  —1  ilarvL^  for  — for  ^oP 

-Urr-ot  for-  ^ r^o  for 
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Integrating  Eqns  (112)  and  (113)  with  gives 

= i.  Cu^e  t CC^ 


(116) 


where  is  the  value  of  at 


and 


=■  ^C-ciL^t  *<ic^j± 

^iCLL^i'^cr  t +c, 


where  C2.|  is  the  value  of  /)l^  at 
solving  for  "t  in  (116) 


(117) 


t = 


substituting  't  into  (117) 


'K.&)  -CC^ 

^ O (Xrvi^ 


(118) 


for 


^tClL^  C'A'i - acc,/)i-.fe)+ cC)  (119) 

+cc.r^e)-<r,)  +C| 


-¥,t)  =_Z — -C, -y,^)  -f'CCr  +C,^A)  -CC^  -f  C, 

^ ^(JryrM^  U/t*'**^  ^ U'n^ny  ^^yv\dLpe 


where 
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(121) 
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for 


t^=-a 
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the  switching  curve 

Is 

- 1 /Xjd:)  1 

S.CCL^ 

(124) 
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o 

(126) 

It  can  be  noted  that  this  controller  design  Is  almost  completely 
Independent  of  the  aircraft  type.'  In  the  low  speed  range  where  aero- 
dynamic controls  are  not  effective,  this  design  will  help  stabilize 
the  roll  mode  of  any  aircraft.  The  only  relationship  between  the 
aircraft  and  controller  is  that  the  thruster  force  is  a function  of 
roll  Inertia  and  wing  span.  Thus,  this  design  becomes  very  versatile 
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and  applicable  to  stabilize  the  roll  mode  of  any  ACLS  aircraft. 

A somewhat  similar  analysis  will  be  made  for  the  controller 
of  the  yaw  thruster.  The  criterion  for  directional  control  is  to 
keep  the  aircraft  on  the  runway  centreline  during  takeoff;  this 
can  be  accomplished  by  minimizing  the  lateral  deviation,  y,  and 
rate  of  deviation  from  the  centreline,  y.  This  deviation  and 
rate  will  be  minimized  by  yawing  the  aircraft  In  a direction  to 
oppose  the  disturbance  with  the  use  of  the  yaw  thruster. 

Prior  to  the  Installation  of  the  air  cushion,  the  direction- 
al stability  of  the  Jlndlvlk  was  controlled  by  a batsman  at  the 
end  of  the  runway.  His  job  was  to  steer  the  dolly  (Fig.  1)  to 
keep  the  aircraft  on  the  centreline.  The  same  batsman  will  visually 
sense  lateral  deviation  and  deviation  rate  and  control  the  yaw 
acceleration  to  Indirectly  control  the  lateral  deviation  and  rate. 

Assuming  that 'the  pitch  and  roll  angles  are  kept  small,  then 
the  lateral  acceleration  (to  correct  a lateral  displacement)  In 
the  Inertial  frame  of  reference  of  the  runway  Is  a function  of 
the  thrust  and  yaw  angle. 


LJ  ^ THPtUST SIN  y 

y (127) 

and  since  V will  be  small  (^30°)  to  keep  the  aircraft  on  the  run- 
way centreline,  then 


^ THPMST 


y 


(128) 


Equation  128  can  be  Implemented  as  shown  In  Fig.  14.  With  these 
control  loops  the  desired  yaw  angle  to  zero  lateral  displacement 
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can  be  determined.  Referring  to  Fig.  14 


the  Inner  loop  open  loop  transfer  function  Is 


^H),  = 


and  the  equivalent  closed  loop  transfer  function  is 

“ I + <siH 
- Kk, 

S’hkki 

the  outer  loop  open  loop  transfer  function  is  therefore 


_ _KXxJ<sl. 

~ SCstXK,) 


(129) 


(130) 


(131) 


(132) 


this  gives  a root  locus  as  shown  in  Fig.  15.  For  a damping  ratio 

of  0.7  the  closed  loop  roots  are  located  at  ^ 5 — which 

gives  a closed  loop  transfer  function  of 


but 
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therefore  equating  denominator  coefficients  gives 


K K,  ==  K-^K. 


or 


f 

'X 


= K, 


since  K|ls  arbitrary,  a value  of  0.3  was  selected^so 


R,  =0-3 


and 


so  the  desired  yaw  angle  Is 

The  yaw  angle  Is  associated  with  the  yaw  thruster  force  by 


Y ^c^Fyr 


where 


C<  2--73  Xid 


in  matrix  form 
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(141) 


This  equation  Is  the  same  as  Eqn.  114  for  the  roll  thruster,  and 
similarly  a bang-bang  control  exists  for  which  the  switching  function 
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Is 


SXYT^  ^ i-  V 


lii 


ac  F, 


/r/iAx 


(142) 


and  the  optimal  control  law  Is 


yr 


-IS  sxyr  >0 

^yrnAx  ^ <0 

-Fyrn/»iC  ^ y >0 

FyrMAjf  , sxyrs^o^  V <0 

o .^y-o^y^o 


(143) 


The  Implementation  of  this  control  law  would  drive  the  yaw  angle 
and  yaw  rate  to  zero  In  the  minimum  time,  but  the  directional 


control  problem  requires  that  the  yaw  angle  be  equal  to  the  de- 
sired angle,  This  can  be  accomplished  by  shifting  the  switch  - 

Ing  curve  by  the  amount  . 

sxYr=  r-Vi 

This  change  In  the  switching  curve  means  that  the  yaw  rate  and 
the  quantity  ( ) will  be  driven  to  zero  In  the  minimum 

time,  or  T will  equal  Yd- 


0,  C 


(144) 
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Chapter  V 

The  Computer  Program  and 
Simulation  Results 

General 

T.*>  EASY  Dynamic  Analysis  Program  to  Aircraft  Modelling  (Ref.  5) 
formed  the  major  portion  of  the  computer  analysis  and  simulation  and 
the  air  cushion  system  was  modelled  by  the  Foster-Miller  model,  as 
described  in  Chapter  III.  The  computer  program  is  listed  in  Appen- 
dix B.  In  brief,  the  EASY  program  provided  the  means  for  an  analysis 
of  the  six  degree  of  freedom  rigid  body  dynamics  of  the  Jlndlvik 
drone  and  the  Foster-Miller  model  estimated  the  ground  forces  and 
moments  transferred  by  the  air  cushion  to  the  airframe.  Additional 
FORTRAN  was  used  to  reflect  the  Jlndlvik  autopilot  and  the  designed 
roll,  pitch  and  yaw  concrollers  in  the  EASY  program.  Simulations 
were  performed  to  obtain  time  history  comparisons  of  the  uncontrolled 
and  controlled  aircraft  models  with  a crcsswind  driving  function  and 
an  initial  pitch  angle  to  simulate  flying  off  of  a 2 inch  step. 

The  Computer  Program 

Figures  16,  17,  and  18  show  the  computer  block  diagrams  of 
the  aircraft  dynamics,  the  longitudinal  autopilot,  and  the  lateral 
autopilot,  respectively.  Understanding  the  symbology  used  In  these 
figures  would  require  considerable  referral  to  Reference  5,  but  a 
general  description  of  the  schematics  will  be  given  here. 

In  Figure  16,  SD  performs  the  six  degree  of  freedom  rigid  body 

dynamics,  AV  calculates  the  aerodynamic  variables,  LO  calculates  the 

longitudinal  force  and  moment  sum,  and  LD  calculates  the  lateral  force 

and  moment  sum.  The  terms  FX3S2,  FZ3S2,  TY3S2,  TY3S2,  TX3S2,  and  TZ3S2, 

shown  feeding  Into  LO  and  LD,  are  the' sums  of  the  engine  and  external 
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Fig,  17,  Block  Diagram  of  Longitudinal  Autopilot 


Fig,  18,  Block  Diagram  of  Lateral  Autopilot 
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(l.e.  air  cushion  system  and  controller)  forces  and  moments. 


I 


The  longitudinal  and  lateral  autopilot  functions  shovm  In  Fig- 
ures 17  and  18  were  developed  from  the  elevator  and  aileron  transfer 
functions  given  In  Reference  3*.  The  maximum  deflection  of  the  con- 
trol surfaces,  the  gearing  ratios  between  control  surfaces  and  ser- 
vos, and  the  maximum  slewing  rates  for  the  servos  were  also  programmed 
into  the  model  by  FORTRAN. 

Air  Cushion  Program 

The  air  cushion  system  was  programmed  with  ten  subroutines  to 
the  EASY  program  and  the  flow  chart  of  the  subroutines  is  shown  in 
Figure  19.  The  functions  of  the  ten  subroutines  are  as  follows: 

FM  is  the  main  subroutine  which  calls  and  interacts  with  the  re- 
maining subroutines:  it  also  determines  the  appropriate  fan  curve 
and  contains  the  integration  routine.  HC,  TK,  SE,  CO,  PR,  CL,  SI, 
and  SP  form  a set  of  subroutines  which  need  the  aircraft  position, 
cushion  and  trunk  pressures  and  ground  profile  as  input  parameters 
and  they  calculate  various  areas  and  volumes  associated  with  those 
parameters.  HC  calculates  the  value  of  side  trunk  height  for  a given 
cushion  to  trunk  pressure  ratio.  Subroutine  TK  takes  that  height 
and  calculates  trunk  cross  section  parameters.  From  these  parameters 
SE  updates  the  trunk  division  parameters.  Subroutine  CO  transforms 
position  vectors  for  each  trunk  centre,  from  the  vehicle  frame  to  the 
ground  frame,  and  then  calculates  the  distance  between  each  of  the 
trunk  segments  and  the  ground,  and  it  also  calculates  the  ground 
coordinates  above  which  each  of  the  segments  lie.  Subroutine  PR 
determines  ground  elevations  (Input  by  user)  corresponding  to  the 

V 
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ground  coordinates  generated  by  CO,  and  subroutine  CL  determines 
the  hard  surface  clearance  for  each  segment  using  (a)  the  ground 
elevation  value  and  (b)  the  distance  of  the  trunk  segment  from 
the  ground.  Finally,  subroutine  SP  calculates  values  of  differ- 
ent areas  and  volumes  for  each  segment  and  adds  them  together 
to  give  total  areas  and  volumes. 

In  addition  to  these  seven  subroutines,  FM  also  calls  sub- 
routine ST.  ST  determines  the  value  of  fan  flow  for  a given 
value  of  fan  pressure  rise  and  also  calculates  the  forces  and 
torques  for  a given  ACTS  orientation.  Subroutine  ST  also  In- 
corporates all  the  system  differential  equations  so  that  the 
value  of  the  state  differentials  can  be  updated  each  time  it 
is  called  by  the  Runge-Kutta  Integration  routine.  The  forces 
and  torques  calculated  by  ST  are  passed  to  the  EASY  program 
via  FM. 

Simulation  Results 

The  results  of  the  simulation  to  test  the  controller  de- 
signs of  Chapter  IV  are  shown  In  Tables  I and  II.  For  initial 
conditions  of  -1°  In  pitch,  0*’  In  roll,  and  a constant  AO  ft/sec 
crossvlnd,  the  uncontrolled  model  induces  a roll  angle  and  pitch 
angle  that  are  lightly  damped  In  comparison  to  the  controlled 
model.  Also,  the  restoring  yaw  angle  Is  less  than  Che  controlled 
yaw  angle  and  the  lateral  deviation,  y.  Is  subsequently  greater 
for  the  uncontrolled  model.  The  l.S  second  simulation  for  the 
uncontrolled  model  required  15,000  seconds  of  computation  time 
and  consequently  the  simulation  was  not  continued  to  the  extent 
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Table  I 

- SlmulaClon 

Results 

Uncontrolled 

Model 

1 

LATERAL  ' 

TIME 

PITCH 

ROLL 

YAW 

ALTITUDE 

DEVIATION 

0.0 

-1.000 

0.000 

0.00 

2.60 

0.00 

0.5 

.305 

.113 

-.63 

1.52 

3.84 

1.0 

.451 

.474 

■1.94 

1.02 

8.78  ^ 

1.5 

.377 

.142 

■2.51 

1.25 

14.87 

i 

Table 

II  - Simulation  Results  Controlled 

4^ 

Model 

TIME 

**  J 

LATERAL  j 

PITCH 

ROLL 

YAW 

ALTITUDE 

DEVIATION 

I 


t 


( 


t 


! 

i 
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Chapter  VI 

Conclualons  and  Recommendations 

The  National  Aeronautics  and  Space  Administration  (NASA) 
has  accepted  the  Foster-Miller  program  as  a valid  air  cushion 
model;  however,  in  the  course  of  the  analysis  of  this  thesis, 
several  problems  arose  which  prevent  this  program  from  be- 
coming an  effective  design  tool.  The  primary  problem  is  the 
excessive  computation  time  required  for  dynamic  simulation 
when  It  is  incorporated  with  the  EASY  Dynamic  Analysis  Program. 
The  Fourth  Order  Runge-Kutta  integration  routine  used  in  the 
air  cushion  model  requires  a time  increment  of  0.001  seconds 
for  numerical  stability.  Since  this  Integration  routine  is  the 
prime  reason  for  the  excessive  computation  time,  it  is  rec- 
ommended that  the  routine  be  changed  or  augmented  to  reduce 
computation  time. 

The  air  cushion  model  assumes  that  the  trunk  is  an 
elliptical  shape  rather  than  the  actual  shape,  in  which  the 
aft  end  is  10%  wider  than  the  fore  end.  This  discrepancy 
impinges  on  trunk  and  cushion  volumes  and  areas,  pitching 
and  rolling  moments,  clearance  and  gap  areas,  etc.  In  other 
words,  it  requires  considerable  evaluation  and  extensive 
modification  and  verification  of  the  program  to  change  trunk 
shapes.  It  is  rec&uariended  that  Foster-Miller  Associates  be 
asked  to  modify  their  model  to  accommodate  different  trunk 
shapes  as  future  designs  may  require. 
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The  air  cushion  model  has  no  provision  to  orient  the  trunk 
orifices  ocher  .than  perpendicular  to  the  trunk  surface.  In 
fact,  the  Jlndlvlk  trunk  orifices  are  drilled  Inward  at  a 45** 
angle  to  produce  more  cushion  pressure  In  the  region  of  trunk 
contact.  Some  adjustment  should  be  made  to  the  model  to  allow 
this  orifice  orientation  as  a design  parameter.  Also,  the 
model  uses  a single  curve  to  describe  the  fan  characteristics 
of  outflow  vs.  drive  pressure,  but  Che  actual  characteristics 
dependent  on  more  variables;  hence,  a fan  "map"  Is  required 
• f replace  the  single  curve  and  adequately  describe  the  fan 
during  all  phases  of  Its  operation. 

A weak  part  of  the  computer  simulation  is  the  evaluation 
of  the  Jlndlvlk  Stability  Derivatives  due  to  the  fact  that  the 
static  wind  tunnel  data  was  extrapolated  from  the  Recovery  Trunk 
Data  and  was  suspect  from  the  beginning  of  the  analysis.  Conse- 
quently, It  Is  recommended  chat  wind  tunnel  tests  be  conducted 
in  a moving  belt  tunne  ith  the  takeoff  trunk  and  with  measure- 
ment of  Che  rate  variables  p,  q,  and  r.  Barring  this  option, 
the  development  of  the  derivatives  should  be  reviewed  and  ammend- 
ed  with  the  use  of  more  sophisticated  data  reduction  techniques. 
Once  the  computer  program,  shown  In  Appendix  B,  is  changed  Co 
encompass  the  previous  recommendations,  it  can  be  used  to  de-- 
fine  the  following  parameters: 

a)  Operational  limits  and  directions  of  crosswinds. 

b)  A "ground  roughness"  criteria  above  which  the  aircraft 

becomes  unstable.  v 
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c)  A flap  deflection  schedule  to  provide  minimum  takeoff 
distance  within  pitch  stability.  The  present  two  flap 
settings  would  provide  a step  Input  to  pitch  and  hence 
should  be  changed. 

d)  All  of  the  above  for  different  vertical  thruster  sizes 
and  locations. 

This  thesis  has  integrated  the  EASY  Dynamic  Analysis  Program 
and  a truncated  version  of  the  Foster-Miller  air  cushion  model 
to  simulate  an  air  cushion  vehicle  during  takeoff.  During  the 
process  of  that  integration  and  simulation,  some  major  deficiencies 
la  the  Foster-Miller  model  ua/e  been  highlighted.  This  thesis 
has  also  developed  and  depionstrated  a technique  to  control  bang- 
bang  thrusters  on  the  wing  tips  and  a bang-bang  thrust  deflector 
on  the  tall  section.  Complete  verification  of  the  controller  de- 
sign was  not  possible  due  to  the  large  computer  resources ' that 
would  have  been  required,  but  the  results  do  show  the  control 
trends  that  are  expected.  The  application  of  wing  tip  and  yaw 
thrusters  to  other  air  qushlon  aircraft  should  provide  comparable 
results.  Also,  these  thrusters  could  be  used  on  Vertical  or  Short 
Takeoff  and  Landing  (V/STOL)  aircraft  because  these  aircraft  also 
have  marginal  stability  and  require  control  enhancement  in  the 
low  speed  range. 
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Fig.  A-1  Lift  Coefficient  Versur  Angle  of  Attack 

Fig.  A-2  Pitching  Moment  Coefficient  Versus  Angle  of  Attack 

Fig.  A-3  Drag  Coefficient  Versus  Angle  of  Attack 

Fig.  A-4  Side  Force  Coefficient  Versus  Sideslip  Angle 

Fig.  A-5  • Side  Force  Coefficient  Versus  Sideslip  Angle  - For  Angle 
of  Attack  » 0.1  Deg 

Fig.  A-6  Roll  Moment  Coefficient  Versus  Sideslip  Angle  - For  Angle 
of  Attack  * 2.2  Deg 

Fig.  A-7  Roll  Moment  Coefficient  Versus  Sideslip  Angle  - For  Angle 
of  Attack  4.3  Deg 

Fig.  A-8  Roll  Moment  Coefficient  Versus  Sideslip  Angle  - For  Angle 
of  Attack  a 6.4  Deg 

Fig.  A-9  Yaw  Moment  Coefficient  Versus  Sideslip  Angle  - For  Angle 
of  Attack  *0.1  Deg 

Fig.  A-10  Yaw  Moment  Coefficient  Versus  Sideslip  Angle  - For  Angle 
of  Attack  a 2.2  Deg 

Fig.  A-11  Yaw  Moment  Coefficient  Versus  Sideslip  Angle  - For  Angle 
of  Attack  a 4.3  Deg 

Fig.  A-12  Yaw  Moment  Coefficient  Versus  Sideslip  Angle  - For  Angle, 
of  Attack  a 6.4  Deg 
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‘ ir(69SI3£).CT.30.l  CO  TO  2000 

CY  — it.  »E-5»3-**3*1.6fE-S*3-*»2-l.  33r-2»BE*2.5&E-3 

CPOLL»-l.  1*«E-3*9E.2.3S5H-1 

_ CO  TO  2000  ; 

CO  CVAW.<..42e-!<*aE-.  00  25 

iFiiesiati.cT.io.i  CO  to  2000  

. CV--I..  38£-5*9£**3*1.66r-3*3i**2-1.13E-2*e£*2.$6E-J 

CR0LL»-1.  Hl£-3*eE-2.3S5E-3 ■ 

CO  TO  2000' 

C 4NCLE  OF  »TT»C<  IS  LESS  T«AM  3.0 

iso  CO*CCl*l-S.li.t-c*49Sie£»**3«3roit-<.*»BS«B-)»*'2-2.'37E-<.*A8Sia£M  *0 

XOSI9P)  

lF(A2S(a£I.C£.60.l  CO  TO  2000 

IF<a£.LT.-».l  C0_TO_50 

IFCJt.Oi.T.l  GO  TO  60 

CYAm-7.0E-3*a£‘*5-1.0T;-6»3E*»<.*l.»5E-S»BE**3»7.<.2E-S»9E**2-1.13E 

*-3*8£-1.3<.>3 

CTsS.  03£-7*3-*»5*2.  lAE-6»e£«*4-9  .0  3e-«*BE»»  3- 1. 6<.E-i.»aE*«2-7.  21E- 

KJ*ec*i.l6£-3 

CPQLL.-l.  5Jl--3»6£-i;26£-3 

CO  to  2000 

_Jfi,  CTAH«.,3i-i*5£*3,67£0 

■"■lF(A2S(3£>.CT.33.t  CO  TO  2000 

_____  CT«3.01-.7*3«:*»5*2.  1I.E-6«B£**4»9  .a3E-5*eE»»3.1.6tE-»»9£»»2»7.21E- 
Xl*e£*8.l4£»3 

COOLL»-1.5  31F.3*BE-1.26£-3 

CO  TO  *200  0 

crAK.L.i'-i'Oi-s.ooiri  _ ■ 

lF<A;S(3il.CT,30.)  CO  TO  2000 

cr.3.  C 3£-  7 ‘O-  ••  3*  2,  l.£.A«B- «».-9  .0  3E-5*BE  * * 3»  1.  6>.E-i.»BE»»  2-7.  21E- 

*■  X3*3E*«.16e-3  *■ 

C »OL  L • - 1 . 5ji;  £-3*BE-1.26£-3 

CO  TO  2000 

_5  ANCLE  OF  ATTACK  IS  LESS  fhaM  5.3 

390  C0*C01.«-<..12£- j*49Sl3£l**3*2.5£-<.»46sle£)**2-2.3LE-L*A9SieEll*COS 

XI8R)  

lF<A5S(3Et.C£.60.)  CO  TO  2000 

1F15-.lt. -7. I CO  tq  70 

IF10E.6-.7.I  CO  TO  00 

CTAH.-i..li.:-7«a-*«5-7.67£-7»BE»*F»«..36E-E»B£»»3*6.5<.E-5»a£»»2-l.<.7 

Xi-3*9E-l.OJX-5 

C7»2.9  3£-6*T£**3>3.72E-7»BE»»i.-2.  03E-<.»aE«  » 3- 1.  9£-5»BE  »»  2 -<..06£-3 

X*e6*l.  15£-2 

CR0LL«-1.5*7E.3»5F-t.  555E-3 ; 

CO  TO  2000 

70  CTAM5i..;£-..qv.4.g2£-3 

IFIAeS (9£ ) .CT.  30.  ) CO  TO  2300 

CT.2.<>3£-t*a6*‘5.3.  72£-7»e£»»i.-2.33£-4»0£».»3-1.9E-5*9£*»2  -<>.06£-3 

X*9£*l.36£-2 

CP0LL»-t.9*7E.3«9:-l.E55E-3 

CO  TC  2000 

_aO  CTAM.<..SE-.*3£-3.05£-3 

IF1A9S  nS)  .CT.  30.  I CO  TO  2030’ 

CT«2.9  3£-6*3  :**i*  3.  72£-7*a£*»4-2.  83£.i.*3C  * ♦ 3- 1. 9E-5*BE»«2  -I..06E-3 

X*8£»l.l6E-2 

CP0LL.-l.977E- 3*  9E-1.555X-3 ; 

CO  TO  2000 

C ANCLE  OF  ATTACK  IS  CaEATE*.  THAN  5. 

■'195  COiCOl  .12.  1S£-5*40S  (OEt  l.l.Cli-L’ABSME)  ••r-1.2'*E-3*A5S  IBEIl'toS 

X ( PR ) _ _ 

IF(A9S(3£t.',E.6(/.t  CO  TO  2800  ' “ ' ‘ 

lF(n£.LI.-7.1  CO  TO  90 


'v■^\  AVA'.1'w'.£  COfV 
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JH9i.Ci.7,l  SO  TO  100 

_ CT»i(.l,?8r-5‘''£‘*3»l.?3E-5«l£**9-6.21£-<.*8t-5.33E-<.  

CT*t.76--6*3£**;-l  .91J-7*9£**<i-1.6E-»*9E*»3*2.7?E-5*B£**7-SV»ZE^5 

**nE*l4  fli*£-z 

CKOLL'-Z.  39ZE-3*9E»i;i.hSi-3  “ 

_ CO  TO  ZOOO  

90  “ c7AM»-.2»£-<.*9E».  o'ous  : : 

lF(»aS(9£t.ST.30.t  CO  Tfl  2000  

C»»l.76i-£*3i»*5-l  .9l£-7*9i**i.-l.S£-«i»3E»»i»2.72E-5»e'E'»»2-i.»2E-3 

• X^PE.l.OxE-Z  

CXOLL»-2.  302E-3*8E-l.»i.9i-3  ' 

; ^CO  T0_2000 ■ 

lOa  C»Au.i..2i,E.i,*9E-i..3£-S  ^ 

_IF  U5S(9£1  .CT.  30.  t CO  TO  2000 

CV«1.7e£-6*9i*»j-l  .9l£-7»3£**<,-1.0£-<.*OE**3»2,72£-5“*eE**2-5.»2E-3 

*»9£M.0i.£-2  „ 

CSOLL«-2.  302£-3*8£-l.*.>.5?-3 

_.i000.  CONTINUE ■ 

C OE^IVATiON  f]A  OIH;Nrio‘i»L  OiBIVATIVES  ISTABILIT  Y AXISI 

2MOS«-THQ»3«C/i.. ‘CLAOR 

20S*-0S/WT*C/2.‘CLQ8  • 

2C£LS*-0S*CL£L9  

«N0S.Et-0*S*C‘*2‘C‘HA0R/<,. 

HOs»is/vT_*:^ Z'CMOA.'i! 

HOiLS«CS*C’1ELR*i0. 

_ YORtOS*  CT  

KCR«OS»9*CYAi4 

tCo«C$n'CROLt 

“ ■ ■■  Y<\S*RHO*S*1/<t. ‘'CYRR  ' ' 

YPS*RhO*S* 3/A. ‘CVOR 
«RS«OS/VT*  3**2/'..  •rsPR. 

KRS»CS/YT«9»*2/A.»CNP£ ; 

LRS«0S/VT«3»*2»CI.»»./2. 

tPS*3S/VT*  3»*2*:i,P^2t 

i.oEts»os*e«ci.AH. 

XOS«0.  $XW3S«0.  IHOELSO. ; 

C derivation  of  OlNfNTIONAL  OiRIVATIVcS  I9C0Y  AXIS  FROM  STA9IL1TT  AXIS) 

XO«CL*OS*StS(AR»-CO*OS*COS(ARI 

20**CL*0j<:33 (AR|-CO*OS*SIn(AR) 

tl0«CN»CS»C 

xwo*  ( XNOE*  :OS  < ATI  *•' Z-ZMOS'SINi'AR  |■•C^S^AR)  ) •COSCBR) 

2A0m|  itoS'OOS  (AR)  ••2.XMDS*SINIARI  •COS(AR))»COSC0R) 

r1A0*NM09*C3S<ARl*C0S<?RI 

YP«YPS*005(ARI-YRS*SIN(ARI 

YR»YeS‘COS  HR  I »yp3*SIN(AR) 

X0E»-20ELS*SINiAP|  •COSORI 

XC«  »X0S*CC3(API-70S*SIN(AR»  i•Co■s  (SR) 

20«  »70S»CO;(1PI •YOS«SIM(AR) I • cost  SRI 

20i«2o:LS*:o5(ARi»:osnRi" 

tP«  (LPS'COKAR)  ••  2-  (l.P3«uPSl»SIN(AR)  ACOSIARI  ♦NRS’SINIARI  •*2)  'COS (8 

XR)  ■ ' - - 

I.OA«(LOtLS‘COSTARI-NOtLS»SIN<AP(  l•COS^ORI  __  

tR»<(.RS»C03(AP|»*2-  (nP3.lpsi*3IN»ARI  •COSIAR)>nPS*$IN(ARI  ••2I»C0SIB 
XR1..  . 

MP>  <NPS«C0SIAR)*'  2*  INP.$-LO$l  •s’IMTArV'C'oSIARI-LR'S'SINIARI  ••21  'COStB 

. XR)  

NR«INRS*C03(AR)»*2* (LRS»NPSI •$ IS < AR) ‘COS f AR) «tPS*S IN (AR) • *2 1 •COS«B 
XR) 

nOxNQS  INOEjhOECS 

_20_  _CONTIMU- 

XO  lO*xo' l<0El0»ii0£  170  LO.20  3>i0L'0'2A0'?  20  LO^IO  T20ElO*20C 

MO  tO'MO  SYAOLO»NAO  IMO  LO*nQ  SICCLO.MOt  SYR  ID*YR  SYO  l.O»Y8 

YR  LO'YP  'tc  LO»l,')  ILP  LO»l.P  JLR  LO"lR  5l.rAi.0«l.  DA  ANS  LO*NB 

NP  LD*NP  hr  lO«N»  103  AV.03  ScLlUO'ClE  IAirSO«Al.T  SAL  AV«AL 

UE  AViO*  <T0  AV^PO  iFC  AVsPO  IVT  AV'VT  iOO  AViOO  INSP'12 

tOR10«lOR  3TOPtO*YOR  IriORLOiNOR 
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LOC*TION«liiZ,iH 
rORTRAN  SIATEHiNTS 
RPH»J  JO— 

KIH 

10  call  p*-ihy  ri.pjM 
_ in.css‘’'',RO'<.  v:s 

i ALTsa.ROwSO.Pi r 
X SO.C  SO.«  so. A 
*T«.L2  TX.PHIIX, 
XFf.OTAFM.ATRCH.AT 
XSHA°ct 

_ 

SX«ROLSO» l»0  AVI  • 

IFISX.GT.O.ITXCON 

IF<Sx.lt. O.ITrCON 
IF(  (SX.EO,  I).  t.  AND 
1F( «S«.iC. 0. I . ANO 
1F( (A35(SO.SOI  .'0 


•TH  rM-l.«-.E-2*TH  EM*»2*o.2JE-6»tM  EN»*  J- 1. 6TE-»» 


F1.PT<p*1.o»(.PH,WTXS2.PFMFH,0FNFM,QFXF1,  IPPF' 

F.i,pH0P1,0VKFH,  VPlPM,OV?rH, OVCF>A.lXtFW. 

SO.PHrOOr.THEOOJ .VCHSE.ACHSJ.U  SO.N  so.al  av.p 
P.MSS.PAt,  TEn.MYI,  OPl.Rl  TX.PZ  TX.Ll 

IX  ST. ft  st.iz  st.tamso.xz  it  j.optfh.otc' 
KAT^.ATKChC.ATKATC.QChAT.V  S3.PM2TX,L1IS1,l2IS1.1 


POLL  COSTPOLLER  ♦«»♦.»>*»«»»««»»«»»» 

APS(’P0  AVI/l2.*L'KAX*0.1tE-3t 

a-UTAX'ZO. 

= UKAX»20. 

. (“o  Av.  r,r.  0. 1 ) txconj-uhax»2o, 

. (“0  AV.  LT.O.  M TXC3H=L'MAX»20. 

.O.Ol.AMS.  lAJSdO  AVt.EC.O.O))  JxcONfOi 


£••>••«•••••••>••••  YA 

UTTIAX^IH  VUii, 
SIECp-O.IO  ‘TO  so- 
STTX=4tVS0-SI£^P 
tF(STIx.C  T.O. I TIC 

if(sttx.l:.o.  i'iic 

1F((STTX..'3.0.).A 

IFIISTTX.tC.O.l.A 
1F( (TAVSO. iO.  3!- 
PiTCH  CONT 

SXPtrCH.PITiOMJO 

■ IF($XFITCH.OT.O.  I 

iF(sxPircH.Lt. 0. I 

IFOSTPtlC-.iO.C. 
_ IP(  (SXFtTC-i.tO.O. 
IF( (iss(Pirsoi.€0 

i 

Txjsjtrxcov  jTTjs 

TT2SS»TJ  ST 

LOCATIOLyIJt.SJ. tM»«IS» 
L0CATI0S»12.10.1 V=JTS»A 
tCCATIOK  = »L.L3.  I (p'JTSsA 
_tOCATIOF«3I.:3. INPUTS»L 
"fOPTRAK  STATCIiNTS 


H COMTpOLLiR. 


XZ  IT  1 

SO*AOSCP  sb)*18S.  I'/LPTTtiAX 

0‘I«-UTTHAX«  10. 

OH=  UYTMAi'lB. 

M3,(f  SC.'.T.O  . ) I I2C0H  = -UYT1AX»10t _ 

M3.<*  SO.lT.O.  ) I T2CON«uYTMAX»10. 
Cl.ArO.  (P  S3.-C.0.  HTICONyO.  _ 

ROLLUP 

AV)*A1SI30  AVI /<2. »UHAX»6.$2E»1) 
TYC0Ha-U'HAX»l2, 

TYCONa  UM4«a»2.  _ . 

» . INC. (QO  AV.CT. b;) JTYCCNa-UPAXPie. 

t.AN0.(O5  JV.LI. C. » tTTCOHa  UMAX'lE,  

. 0.  01 . AaiC.  <A3SI0a  AVI.EC.O.OI)  TYCONiO. 


JaTYCO'4  ST2J>3aIZCON  lF22S3a-FY  ST  STX2S1 

V.SJ.MC  t(XaELc) 

V.LO.SJ.HC  ZIXaAtLI 

O.LO 


iLECa.Zja  O SO-0. I a.SJJ*  tPITS0-0.» 

_ XI  IFaJccc : 

tOCATIOXalAP.TF 

LOCATIOAalL7.:i  r.ISPUTSaTF 

'LOCATIOMieT.'';  1 , iNOiJTSaTF.ir  7. IT  » 

FORTRAN  STATEY^ntS  - 

IFIXy  hc  1.3E,10.I  XL  lie  1 = 10. 

IFIYL  pc  l.LE.-lS.t  XL  YC  la-15. 

tOCATIOAai;?.!!  l.iNPUTSaMC  I 

.LOCATIOAaSI.IT  J.INPjTSaSOfXft.XJ 

L0CATI0K=39  .II  l.lHPUTSaSOIYOaX) 

FORTRAN  STATIpenTS 

AILCa  .19EMP  S0-0.>..;.»*(POLS3-0.>a;zaR'  so- 
XI  TF  ZaAlLC 
lOCATIONaZLa.fF  2 

tOCATIOLaZLS.IT  J.lPP'JTSaTF  2 ^ 

LOCATIOLa266.MO  I.  IN»'JI  3 = '.f'  ?.'lf  9.  ITIO’ 

FORIRAN  -JIATEY-NIS  

iriYL  -C  Z.Ot'.a.  I XL  PC  2a9i 
IFIY-  YO  2.LE.-A.  » XL  PC  Za-o. 
lOC(.TIOLa25L.m  I.iNP'jTSaPC  2 
CNO  OF  HOCEL 


BSi 


COPY 


i 

i 
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rOKTRAN  STATCItHTS 

SUiiROUTINr  rHMT,ocs,PT<,e3LtVT<,PCM,0»’N.or*,  IPP.CSI.ePttVCStR 
tM0,0V<,tf»t.1VP,Dvr,S<T,*C6.rH:,OM<-,i>H£00t,TM£00T,vCHSZ,4CMS2  .U 
• i.M.iL •P.o °AT.  ttH.HY:.  OPi.Bi  tk.r?  tk.li  t«,i 

K2  TK.PMITH,  TX  ST.FY  ST.T7  ? 1 ,S1E  . ICC . OPTFM  .OTCPN.  aTAFfl , ATCS7.' 

XATASJ.ACCS  J.AAC1?,'1''.‘<«  1IS1.(.?:C1  .ISMAPEI 

REAl.  l.LS.  HASS.  (.I.ttl*  l.LlT<,L2:Sl.ieT<,HSS 

OllliNSlON  TCIJI.SI  (t  JI.O'LTltl?  ).SL<.C0(27  ».2C*SEIJ2  ),rriS£ISi  I, 
K.OFLSCI5?  »,ITrSE»32  l.2SIS£<3;  >.lSGSe<3Z  *,VTIS1C3Z  I 

X.ACISKSZ  l.YGHCt(S2  1.ACPS2H2  ».AIIS2I3?  >,PERS2I1Z  I.ATgSZ(3Z  » 

t,XCMS2(!2  l.*T<S2M2  » . 2CHS2  < 32  I . 3Tk£2  I 32  3.01  ST<13).TC  PRC52  ). 

XTO »133 ,11 1 131, r;( 1 33 

data  0I'.RA3IAn73.  l<ri592%S3.0.017i.S32/ 

_ CTIr-»,0001  JTlMC.t,  

C AlP  ruSHIOli  LANSING  SYSTEM  ePOGPAM 

HSS»77.»  fl’t  /COI*!.-.  :CAF=.5  *NSP*32  *vco*o.  _scc^i.ir_ 

CG*  1.5  IFF.c.  S®4T*21iS.3  1LS«L.12E  S0=.A17  tA>.6<)S  SC 

£«.52*.SLi.5.  I'  SHYIM.O  :mjl9S  SAm.3.hE-4  S5h«.0>33  SLI»0.«fl 

CGP«.S'  S&£C».l  SC*X*  1.17  JT£H*7b. 

TM£OOT»P  tPMiOOr.Q  SSXT.-n  _ 

HASS>*1SS  ' $TEHOAT«r£H 

^THyAS*TM£  |OmI^»ph£  SO’Vi.ohEOOT  SOTHE TA* THC0(U 

NSTb“»NS®‘  JPp'LNsPPL 

.c  OAT*  ACOUlSfJI’'^ 

VPL««.313  3VFAN«.<.6a 

VPL»W®LN» VF*N  

00  100  l»lVNSfOP 

otLtAdiso.j ; 

too  CONTINUE 

AM0tl.2i.i/H.5Q.0*T£HP4TI 

~HY:HYi 

C HINIPON  trunk  MilCHT  IS  PARKING  5LAT3ER  HEIGHT . 

HY*AHAIl(Mri,70> 

CALL  T<(ISNT<.Ph;T»<.»2  TK.P.l  TK.PHITK.LI  TK.L2  TK,  MY  . A . £ . L .LSI 

ISHAPP.ISMTK 

in  ISHA»E.  iO.  0)  CO  TC  1?9 

lCLEK»d 

_ call  SKITYSi.SI  SE.OI  SE  . 7CXSE.  ICXSE . PELSE  . XCISE.  7C ISE  . ICLFM.R2  T 
iK.PM2Ti<,02tSl.LS.M,N,0,tSCScl 

CALL  SKOZISl.At  Si,  42  31, SI  SI. Cl  Si . OXIS 1 . ? : IS  1 , LIIS 1 ,L  21 S 1 ,RIIS_ 

ll.P2ISl.SrlISl,»H2TK.P2  TX.OI  SE.EE  SE.MY.Rl  T <,  PHI  T< , ISLE  H ,L  1 Tk.L 

X2  TX.I lYs: ,0, A ,E ,N, 1, VT ISt .4CISI .LSI 

pPLH»i7.  0 23-3.  l7:-2*PF'1,i.;3£-6‘P  PH»»2-2.  d*G£- 1 0»RPM»»  3 

_ PPANiPPLH  _ 

C INPUT  PRESSURE  ANO  OUTPUT  FLOM 

in  (PPM.CE.  119  17.  ) , AND.  I’PN.LT,  12075,  n C0TO»7  

IFI  (PPR.Of  . 12075.  I , AND.  IRPi.LT.12715.  1 1 CO  TO  A* 

in  ISPp.Gl.  11211.  I , ANY.  IRO  i,tY.  12351. 32^  CO  TO  09 

1F<  IRPN.CE  . 1 2351.  3 . AND.  <RP1.LT,i2i.Y9. 1 1 ■■C0''l0'90 

irncPM. G-:. 12-59.  3, AND.  IRFn.lI. 12127.3  3 CO  TO  91 

■ in  (ROM.&e.l2«i27.  3.  AND.  (RP'I.lI.  12765,3  3 GO  TO  92 

_ IF<  loPN.CE.  127G5.  3 . AND.  C^o^.LT.lSOOJ.  3 3 CO  TO  93 

IF  ( IPPN.OE.  12903.  3. 4NC.  <R“-3.LT.120Lt.3  3 CO  TO  93. 

IFl  (OPN.CE.  1 lOct.  3 . ANO.  (RFN.  LT.  13179.  3 )_  CO  TO  95 

in  IR®H. Cl. 13179.  3.  and.  (Rpn.lt.  1 33  17, 3 3 GO  1(196 

_ _ IF ( (OPN.GE.  1 T317.  3 . ANO.  (RPH.LT.I 3t55. 3 3 GO  TORT 

IF  ( (ROM,  C£.  1 iCSi.  3 . ANt-.  (kPY.L'.  I 3‘-93.  I 3 CO  10  90 

IF  ( (Oon.cF  . I 3591.  3 . ANO.  (RF.i.LT.l  3233 , 3 3 CO  TO  99 

IF ( (RPn.CE.I  3731.  3 . AND. (RPR.Lt.l,  000,3  3 GO  TO  101 

IF(  RPM.r.T.lYOOO.  3 CO  TO_'. ^ 

1F(  ROM. LT. 11173. 3 '•.0  TP  2 

IF  ( (APN.cr.  1 1 1 7 3.  3 . AND.  (RPR.LT.  1 1 385.  3 3 CO  TO  82 

IF  ( (RPM.  C£.11185.  3 . ANf.  (RPt.Ll.  1 1E2(,3  3 CO  TO  83 

IF  ( (30N.CE.  HE  2 1.  3 . SNF,  ("PY.l.T.l  3E'>1.  3 3 CO  TO  Ik  . 

IF  ( (RON,  GE  . 1 ICfcl.  3 . 13Jr,  (RPR.LT.i  1 79®.)  3 CO  TO  85 

IF)  )RPN.Ci,ll79l.  3.  AND.  (yp-I.LT.I  1 937.)  3 CO  10  86 
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t PRINT  ! 

3 FOR>'»It*OP^  LON'R  JNO  OF  fANji*®*!-  

CO  TO  199 

I,  PRINT  9 

5 FOSNATI*  off  uP^fr  r>4[)  qF  FAN  MAP»» 

CO  TO  199  

"CFAN  CU»V£  for  ifATIC  ITfRATIONS 

QFANtlS.Oe  -.If!  •»f AN»f.34f. J *PFAN»»g-l,eSE-9  *PFAN 

*••3  ICO  TO  no 

es  QFAN«15.oe  -.3J9 • = FAn»3.#i.E-S  _ «PFAN»»I-1.39E-S  »PFAN 

*••3  ICO  TO  no 

«h QFANMI.  j3 lU? »PFAN»1.7ht-3  «PFAN*»t-6.3e»6  «PFAH 

x»»3  jco  TO  no 

_ 45  eFAN«10_.«2  -.1?^  •^PFAM»1.29EtJ «PFAN*«Z-<..SC£-t  »PFAN 

■*••3  ■ 3Gq'  to  no 

46  0F4N»9.73  -.077 •AFAM»  7.  I5£t  A «PFAN**2-7.6£-6  »PFAN 

*••3  SCO  TO  no 

47_  0FAN*9.97 -.07_7 »»FAN.r.  95E-A  »PFAN»«g-7.A£-6  »PFAN 

'X“3  ‘ iCO  TO  ilO 

_66  QF4N*1£.39  -.iS ^faFAN*  U 32E- 3 «PFAN*»2-3.96f-6  «PFAN 

■ X»«3  SCO'TO  110 

_ 49  _OFAritie.38  _-.  333 ^»F_AN»  2,  S7E^3 »PFAN»»  2-7. 79E-6  *PFAN 

X»«3  SCO’TO  ll'O 

90 OFAN«n,03  -.211 »JFAN*1.7aE-3  »PFAN»»2-a.67E-6  »PFAN 

x**3  ' $C3  TO  no 

„9l  0FAN112.73  -.  1 j •^^FAN*1.0  3E-3 «PFAN»»2-2. 7»£-6  »PFAN 

■■■x-»i  ■ iC3  TO  no 

_92  OF4N«13.70  -.  133 ♦»FAN«  1, 1 7E^  3_  »PFAN»»2-3.01E-6  »PFAN 

' " x»«3  SCO  TO  no  ■' 

JS_  0F4n«13.7A  ’♦IA_2 »»FAN«1.  05E-3  ♦PFAN»»2-2.6£-6  FprAN 

■x»»3  SG3  70  no' 

JA  0F4N«12.  05  -.  006 »t’FAN»i  .SEE- A *PFAN»»2-l.S£-6  »PFAN 

'x7»3  SCO  "TO  no  ' 

9S__  CFAN«12.97  _ -.lOA »gFAN«t.93E-A  *PFAN»»2-1.63£-6  »PFAN 

it»»3  ■■  303  To'nO' 

^0^***«lt.59  _ ♦>FAN>3.7e£-A  »PFAN»*2-9.23;-7  ♦PFAN 

“ X773  ■ 'SGO  fo  I'lO' 

^97_  OFAN.12.6  -.  079  *’»FAN«a.79£-A »PFAN7»2- 1. 09£- 6 »PFAW 

■x»»3  3G3  TO  110" 

_90__  OF4n«10,29  -.2^ »PFAN^1.29E-3 »PFAN7»2-2.6IE-6  FPFAN 

x*73  sgi'to  no  ' 

99 0F4N!l9_i.0  3 -,2\i »7F4W.l.  29E-3  *°FAN»»2-2. 53£- 6 »PFAN 

x»73'~  sc'i  ro"'no 

_10l_  0FAN.32.07  .- . A 2.7 •AFAN.2.9 3 ?.PF AN r 3 . 6 3E-6  *PFAN 

X»»3 

_c  ch4ng:o  0 fpqn  mass  f^om  to  volume  flop 

il'o  QFAN.'QFAN/(RM0732.  2'» 

OF  A NX  r OF  AN 

qfx.o'fanx 

OFN.QFAN ; 

pfn.pfan' 

_ PPL.PPIN  _ _ 

'c'  CONVERT  OEGRFFS  fo' RAOIANS  ' 

TMETAi-TNiTA^-RAD  JAN 

PMIt'iPHlE’TAOiAN 

_ OTH; TAsOTHiTA^RAOIAN  

OPHl.aPMl'TAOIAV  SSIe.';:£»9AOI4N 

ICIN.O  

ICN.IClN 

5 OnTAIN  initial  VAL0f_  OF_OCC'7F  ant  ImITIALIIf  geohetrt 

PCMS5= (=CM. (PTX-PCH I *0. U /7lk 

CSS.PCH9S  _ 

CALL  HCIZ  HG.CSSI  * 

mv«htI*7  mC  

C minimum  tpuNX  N'IGMT  IS  •’A-TSItlG  9LAITFR  HEIcnT 

MV J ANAXt ( mV , , 7AI 
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C»IL  T«(ISmT<,phzt<,B?  T<,*1  Tk.PhiTk.LI  Tk.L?  Tk,HT,*,£, L.LSt 

lSH*B;*ISHt<  

ir(ISM*»£.iQ.0ICO'Ta  199 
ICLFM.l 

■ CALL  SillTrse.BE  SC.O*  S* . ZC X J£, XCxSC.OELSE . *CiSt , £C IS; , ICLFH.RJ  T 

lK.PMZTx,02ISli_LS.>i,H,0.ISr.i£)_  _ 

■ C*Ll  CC«S(.«CO,XCXS;.2CXSt,»r.G,XCC.SIE.PM£.IHE.CC,FF,CGl 

C*ll  PStTG  *91  _ _ 

"C*Ll  CL  IYCHCL',ICN.“HC,TH;,JL<.Cd,XC  PRI 

RhO'’««RHO  _ 

CELL  S1(0?1$1|A1  SI. A?  Sl.s:  Sl^OI  Sl.OXlSl.eEISl.LlISl.LZISl.RlIS 
11.R2I51.SNIS1.P-<2TK.b?  t<.3X  S£.C£  SE.sr.pt  TK.PmUK.ICLFM.lI  TK.L 
tz  TK.  ITt'sE  .O.A.;  .N,  S,  vrts  l.icisi  .LS> 
call  SPtv_T<S2  .WCSS2  .AChS?  .PmPTK.AP  Tn.OZISl.eE  SE.Al  S1.A2  S1.DX_S_ 
It.s:  SI.  ' IirSE.rr.HCL.sr  .tl  S l.  Cx  IS  l . BilSl . XCISE . ZCl  S£  . L2  I 
iX.LlISl.LZrSl.RtlSl  .RJISl.DiLS- .'•JISl.xcxSE.O.LS.AH.MX.MH.Hrl.LX.S 
Sm.VCC.PI  T<.AGos?,4  TAS?,Ar;s2.AC».SJ.  AtCS2.ACCS2.  ISGSE.VTISl.ACISl. 

J(A CR S 2 . A r I S 2j  A T 9S 2 .P£  9E2  . xCh52  . xj 7CmS2. 2TK52  .L  . L 1 TO  

VCMSS«VCHS2 

_VCS*WCHSS. 

CALL  MC(2  *HC.  IPCH/PTXt  I 


HTsHrl»2  MC_  _ 

C NINIRUH  TRUNlt  .<£:CHT  if ‘PARKfuC  UATTER  HEIGHT 

Hr=  AMAXU  Hr,  . 781  

CALL  TxriSHT<<,eH2Ti<,e2  tk.ri  IK.PhiTk.LI  t<.L2  Tk  .Hr  . A.  £ . L . LS) 

lSHAoc.;sHr<  ^ 

IFCSHA®;.  n.o'lGO  TO  199 

call  Si«ITrS;,9£  S-.OX  S- . 2CXSe. XCXSE. D£L5E . XCISE , ZC tSE . ICLFH. 92  T 

tK.PH2T«,D2:Sl.LS.H.M,0,ISGS£l  ' 

CALL  COiSL .CO.XCXS* . 2CX5~ . rCG .XQC.S: E.PHE . TM£ .CC . Ff .CG) 

CALL'PRCrC'PRI 

CALL  CLirCHCL.ICM.PHE.TH'.SLfcCO.rC  «>P>  

call  SKOZISt.At  S1.A2  St.St  Sl.T!  S 1 .OX  ISi  . SEISi  . L US  1 .L  21  SI  • 91  IS 

ll.R2ISl.Sti:Sl.PH2’<.B2  T^.OX  S£.61  SE.HT.Pl  T<  ,Pm1  TK , ICLFH  .Ll  TK,L 

X2  TK.ITrsE.D.A.E.M.H.VTISl.iCtSl.LS) 

call  SP(.VT<S2.VCMS2.aCMS2.»h2L<.F2  TX,021S1jB£  SE.Al  S1.A2  SI. OX  S 

lE.S!  SI.  ITrSE  .rOHCLTMr  .n  S 1 .OXIS 1 . RilSl . XC !5E . 2C is£  . L2  T 

2A.LUSl.L2tSl  ,RUSl.92!Sl  .O-LSi.SMSl.  XCXSE.O.LS.AH.NX.MH.Hrl  .LX,S_ 

3H.VC0.9l  r<. ACPI 2. 1 •AI2.ATC32.A:nS2.AiCS2.ACCS2.:5GS£.VTISl.AClSl. 

XACPS2.ATISI.lTRS2.BE9r2.XCHS2.«TKS2.  ZCmS2.,ZTi<S2.L.LA_T'<I ■ 

OVCHB«  (VCM32-VCHSSI  / ( ('>CM/?TKI*BCHSSt 

OWPsOVCHP 

C»»»**&rNAHTC  si'xijLATrOH* 

INUM*0 ^ 

‘ ovtk^o.oT 

OVKxOVTX  

0VCh.SxT»A0MS2 

ovc*ovcM 

CALL  STtFr  ST.TZ  S T , T x ' ST  , Or.  sY  , f=P  . FP'i  . VPL  . V T XSz',  O' X , 6vH  . P TK  , OVC  . 

trCC.0VP.PC-l._VCHS2.SKT,PH£.OMi,00T.I“;.THrC0T.CK<.PAT,GG.MrI,G£C.HSS 

2.  fiPM,  ■■  ■ - ■ y,  sir  .ohO,AOP52,AChS2,ATAS2.ATCS2,ACNS” 

SZ.AACSZ.ACOSZ.CAF.Cfx.CGO.  ACf £2. ATIS2.ATRS2.PERS2.XCHS2,2CHS2. 

3xTKS?.ziKS».4Cisi.rcMCL.i®TFM.orrFH.cTAFM,ocHtn 

.calculate  ovcM.ovc"P.ovr<  as  cvc.cvp.dvk  .».»».»».»»»»«»»»»»« 
IF(Ini)MI2Q3. 1,200 

_ 1 VCMS*VCHS2  SVT«S»VTKS2  _ _ 

200  call  TK(ISH’K.PH2T<,92  TK.Rt  TK.FhiTK.LI  TK.LE  TK, HY .A.E . L .L$) 

ISHAPE.ISMTK 

lCLfP«l 

CALL  SEIITrS'E.BC  SE,Ox  ST  , ZOXSE,  XCxSE  ,0£  LSE  jXCISt  lZCISJj  lCLiMtR2_T_ 

tK.PH7TK,02ISI  .LS,H.N.t.ISr,S£» 

call  CC(SL  ,00.  XCiS-;  .7CXS£.rcc,xcr  .Sir,PHE,TH£,CC,FF,CC»_  

call  POIYC  pri 

...  call  CLlrCHOL.ION.PHF  .THE.SLLCO.rC  PftI 

f hot  “ ^ T..(j 


call 

r 


lIDZtSl.Al 


Sl.AZ  St. SI  Sl.n  SI.DXlSl.ectSl.LUSl.LZlSl.RUS 
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lltR?TSl.3NtSl.P‘<2T<,P?  TK.'il  SE , P?  ?£,hy,R1  T<,PhiT<,ICLEM,|.1  TK.L 

_ T<,ITrSC.T,».-:.N,'<,v’I<il.iCISt.Lil  

t*LL  SP(VT<S^.v:MC?,^C-C^l»•1^Tl<.F;  rr.0?i;i.9i  SE.il  SI. A?  Sl.OX  s 

It, St  SI,  ItTSi .Yf.HCt  .H,  .ri  Sl.Oxltl.ScISl.lClSE.ZCISt.LZ  t 

?«;,tlISl.t?ISl.«lISl.P?tSl  ,0£LSE,rNISl.«CXSE.C,LS,AM,NX,NH,MfI,LX,S 

3H,VC0,91  TK.AC.as?.  tTAr?,tr:£^.t;,.c?,.^es;,ACCS?.15CSE.VTISl.ACISl. 

»ACRS2.ATISI.AiRS?,PERt2,«CH32,XTKSZ,7CHS2.2TXS2,L,tl  T<t 

CVC*(VCHS2-VChSI/. Ot  ; 

OVCmo  CitC’JtAf  lOM  ••••••••••••••*••••••»•»»• 

VCHSStVCMS?  S=R1T«  ( aCH,«PT<-FCHl  •,i)/FTit 

■ CALL  HCIMX.ORATl 

MYPMX'HYl I 

C HXNIMUH’tRUNit  MtlCHt  IS  PARKIMC  iLAITER  HEIGHT 

MYiAMAXl  < HY,,  78 1 _ 

CALL  T<tlSHT<.PH2t<,e‘7  Tk.RI  TK.PMITk.LI  T<,L2  TK.HY  , A,  £,  L ,LS» 
ISmAPp,IShT<  

iniSHAae.pa.oi  go  to  199 

CALL  St(ITYSi,8£  St,C«  S£ .7CXS£. »CtSL . 0- LSE . YCISE . 2C1SE . ICLPM, R?  T 

lit.PH2T<.0  2ISl.L5.H.H.O.IS&Si) 

_ call  CCtSL  .CO.YTXSE.ZCXSE  .YCG.*CG,$It,PH£,TME,CC,FF,M» 

CALL  P9 (YC  PRI 

call  CL<YCHCL.ICN.PHE,XHE.SL<iCO,YC  PR) 

RhOFH*SHO 

CALL  Sl(02ISt.Al  S1_.A2  St.Sl  SI. FI  Sl.OXIS  1 . PEISl  ,L1  IS  » . L 2 IS  1 . Rl  IS 

ii.R?isi.s*.:si.pH2T<.;2  Tt.oi  s'.re  se.my.p.i  tk.pmitx.’iclph.li  ik.l 

_ *2  K.IlYSC.a.l.-Z.H.H.VTIjl.ACISl.LS)  _ 

call  SF(VT<S:.V'.HS2.tC-iS2.»H2T<.P2  TK,CI21Sl,Bi  Sc.Al  S 1 , A2  SlVO*  S 

It, SI  SI.  IIYSE.xGmJL.hY  .01  S1,0<ISi,?:ISi,xc:s:,2CI5S,L2  t 

2K,LlISl,L2ISl,RnSl  .R2ISI  •3iLSi.tt,lu.xcxSE,0.LS.AH,NX,NH,MYl,LX,S 

3H.VfO,Al  T<,tlP;».ATA'?,  AT:32.AS);S2,AACS2,ACCS2.ISGS£.VTISl_,mili_ 

XACRS2.AT:S’.4r9-:7.p:pr:,«CH52,XIXS2.2CMS2,2TXS2,l,Ll  TKt 

__  OVP*  tVCHSS -WChS?)  / < <FCHFO''<»-»»AT) ; 

PRAT«A-AX1(0.  O.iMItll  (1.  3.  |PCH/PIxY)“) 

(*•••«,,  CALCJtATi  OVT<  ,,»>»»,,«»»«»»♦»»»»»»»»,»»,> 

call  HCIMX.oRAir  

hy.hx*hyi  

e HINIHUH  TRUNX  MlIGHr'lS  parxihs  blatter  height 

__  MY*AYAXl<Mr.*,  78)  

call  T<<IS-*7x,OH2I<rR2  Tx.RT  fxTPHlTx.ll  IX, L2  Tx.HY  ,a,E,L  ,LS) 

XSHAP'tlSHTX  _ 

' IF  iiSHAot.  :a.  ordb  to  hr 

CALL  S3(ITY3r,B£  Si, Ox  St . lOXSE . XCXSf .0£LS£ . XCIGE . 2C ISS . ICLFM.R2  T 

TX,PH2Ti(,a2tSl  .L>.H,H,0.lSGSil 

CALL  CC(SL,C3,Xi;<S£j2CXS£  .Y^G.XCC.SUlPhE.THE.CC.FF.QC) 

CALL  PRfYG  PR) 

call  CKYOhol.ICH.PhE.THE.SLACO.YG  PR)  ' 

KM0PH,RM0 

CALL  Sl(02ISl,AS  SI. A2  SI. SI  SS.fl  S 1, Ox  IS  1 . 9£ 1ST . L J ISl , L 2! Slj RIIS 

11.R2IS1,SNISI.ph2TK.R2  Tx.Ox  S£.£i  St.HY.Fi  Tx . PHI TK, ISLFH , LI  IX, L 

X2  Tx,iTYS£.o.A,:,N,H,v*:si,ic:st.LSi  _ 

CALL  SP(VT<S2,VSHS2.ArH52.>H?Tx.F2  TX.02IS1,9£  SE.Al  si,A2  Sl.O'x  S 

lx, SI  SI,  , IIyS-.yCmCl.ht  ,ri  Sl,0XISl,0£ISl,xciSi,7ClSE,L2  _T 

2X,L1IS1,L2IS1 ,«l ISl  .P?IS1 , 0£LSS, 'MISI.XCXSE .O.LS.AM.Nx.nh.myi.lx.S 

3m,VC0,R1  TY.8G0S2,ATAS7.  AT:s7,ACt.52.  tACS2.AGCS2.ISG5£.VTISl.ACIST. 

XACRS2,ATIS’,ATRS?.°S-S2,XCHS2,XfX52, 2CHS2.2TKS2,L.LT  TX) 

OVX» «V7XS2-VTXSI/.fll 

rOUITH  OROfR  PUH&t-XUTTA  IHTSCRATIOH  ROUTINE 

TI1)»PPL  |Y)2),PCH  $Y»3)"PIX  _ _ 

TIHEifl, 

iP M«0TIM£/2, ; 

00  10  I>1,8 

$TII)*YIII  

YOlll-CT  SKI) 

10  Y (I) ,T(I) ,H«Or  Sill)  _ _ 

TIMPiIIHi »H 

TIlliPPL  IYI2)«PGH  SKDaPTX 


■^51  AV 


" ' I'  A r 


co?r 


82 


1 


I 
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C*l.L  ST(FY  ST.TZ  ST.T»  ;T,D»  <;T.1FP.PPL.VPL.VT'CS2.QF*,OV<.PT<,OVC, 

lTCC.5VP,PCH,V':H5?.S<t,PM-:,»MEnOt  .-H*  ,TM£DOT,C«.P»t.CO.My  I.GEC.HSS 

z,  pPM,  u.  "ir  .PH0,*&PS£.»CHS2.»its;,ATcs£,ACNS 

Si,4ACS2.Ac;52.CAF,CF«,C';P,  ACP'S.ATIi;  ,AI«S?,P£R?2,<CrtS2,ZCHSij 

3x1KS2tZTi<S2iACISl  .t&MCL  .QPT  fh.QT  C F«  ,Ct  AFn  ,0CH4T  » 

_D0  20  1*1, « 

YllI) =OY  ST(I» 

20  ttII-.SttII*M*OY  suit  

Ttll«PPL  »T(21«PCH  IY«3t*PtK 

_ call  STiFy  ST.T2  ST,T»  STOY  ST , 2PP.  pOL  . VPL  , VT<S2.  OFX.  0V<  . PK.OVC. 

lYCC.OVF,PCt,VCHS2  .SST  ,!»H£  ,9H£^OT  .Th-  ,TM£COT,CKK.PAt,CC.HYI,CEC,MSS 

2.  RP-,  L. Stf  ,°hC.ASPS2.ACHS2.ATAS2.ATCS2.ACMS 

32,A4CS2.‘AC2S2.crF  .CFX.CSP.  ~ 4Sf  S2,  AT  IS2.  AtRS2.PERS2,XCHS2 ,2ChS2. 

3xTKS2,2t<S2,ACI£l,15HCL,QPTFM.QTCF»*,OTAF»1,QCiUTJ 

CO  30  I-l.^ 

YZai'OT  STCtl  

30  Y Ct  iSYC ) •3tI1£»0Y  STU>  A 

Y(1I»®PL  JY(2»iPCH_.£Y(3l_t“rK ; 

call  STIFy  ST.Tt  St  .'ty  'st.OY  ST.!PP,PPL.VPL.YTXS2,OFX,OV<.PTK.OVCt 

lYCC.CVO.OCH.  YCMT2_iSXt  .ph£  .OMJOOI  ,THr  .thECOT  ,C<<,PAT,  GG.my  I.OtC.MSS 

2.  ROH,  U,  SIl  .RH0,A&PS2.4CMS2.4TAS2,AtCS2.ACNS 

32,A4CS2,AC;S2,CAF,CCV,CSP,  ASF  S2 , AT  IS2  , A TRS2  .P£  RS2 , XCMS2 , 2CHS2j 

3XTKS2. 2TX5  2.4CI$1  ,YGHCL.QPTF(4,QTCFh,cTAFH.0CM4I) 

: 

HlH/3> 

00  GO  lYl.G ^ 

PPT1»2.0*  (Y1(:»»Y2(IM 

PRT2*Y0(n*5T  bTFU ' J 

~ T<ii«SY<ii*H‘ppri»M»PRr2 

«0  CONTINUE 

fPL*T<11  !“JH*T<2»  3Pt«<»Y(S» 

call  ST(FT  ST.TZ  ST.TX  3T,3Y  ST. JPP, PPL . VPL . VTXS2 . OF x . OVX ,PIX ,0VCt 

tTCG,CVF.PCH,YCMS2 ,S<’  .»HE .pmeCOT , , t m£0CT ,C<< .PA T , GG , MY J , CEC.HSS 

2«  RPM,  I'. Sir  .’mO,AGPS2,ASmS2.ATAS2.ATCS2.ACN5 

32,A4CS2.AC:32,C4F .SFx  .CS®,  ASF S2 , AT  IS  2 . A ISS2 . PE RS2 . xCmS 2 , 2CK$2, 

^5XT<$2.n<S2.ACISl  .YGhCL.OPTfm.OTCF^.CTAFM.OCHAT)  

IF((PPL«ST(in.LT.l,»  GO  To"  60 

lF<Tr~E.LT.flNCl  GO  TO  50 

60  CONTINUE  " 

1NUM«1  

C*»»»*»**CAL'CULAT£  OVCh.OVCmP.OVIX  as  cVc.ovp.ovk 

IF«:NU“120t,ll.201  

11  VCMS«VCMS2  3VT<S'»VTi<S2 

_201  call  TK(ISHTx.CM2TH,a2  TX.Pi  TK.FH1TK.L1  Tk^LZ  TK. HY . A . L ^LSI 

ISM4?6*ISmTX 

ICLFP.l 

CALL  SitlTrSE.OE  SE.PX  SE  . 2CXSE  . XCXSE , DELSE  . XCI  S'F. . JC  I'SE  , ICLFM.R2  T 

lX.PM2T<.02tSl.LS.M,J),0.ISG5- l_ 

call  CC(SL«CO.xcxsE.2CxSE.y;C.x:c.s:£',PmE,‘TME.CC.FF.CG) 

CALL  PPIYC  PRt  

CALL  CHYCmSL.  ISN.PME.TmE.SLI-CO.YG'  PR) 

RmOfmiShO  _ 

call  St(02ISl.Al  "si  .A2  si.si'si’.nr  SI.OXISI.OEISI.ViIS'i.LEISI.RIIS 

^ll  .fi2IS1.3NI3l,PM2TX.Pr  Tk.Ox  SE.PE  SE.hy.Pi  T<  . PMIT  X . iclfm  , L 1 tk.l 

X2  TK.ITYSf.O.A.E.N.M.vTtTi.iClSl.LSI 

call  SF|VT<S2.v:MS2,ACMS2,°H2TR,r7  Tk.02ISI.BE  SE.Al  S1.A2  SI. OX  S 

IL.SI  SI.  ITYSE.y&HSL.ht  ,i:  Sl.CXlSl.DilSl,xcISE.TCIS£.L2  T 

_ 2K.LtISl.L2ISl.RnSl.R71Sl.OEL3E.£MISl.»CXSt.O.LS.AH.NX.NH.MTl.LX.S 

iM.VCO.Rl  T<,AC»S2.1iAS:.lTSSr.A;NS2,AACS2.ACCS2.:s&SE.YtISl.ACISl. 

^XfcCPS2.ATIS2.ATPSr.  P:PS2.Kr.‘i52,x;KS2.7CMS2.2TKS2.l.Li  TK) 

CVC- <VCmS2-VCh3I/.  0 1 

(*••••..  ovchp  CALCUH’TOh  ••••••••••••••••••••••••♦• 

VCmSSiVCkE?  5®RAI » (PCH* FPTK-PCm)  • , 1) /PTK 

CALL  MC(MX.PPAT|  _ . _ _ 

Mr«HX»MYI 

CALL  TK(IShT<.Pm2T<.9?  T<,R1  TK.FMITk.LI  TX.L2  TK. MY . A. S , L . LS» 

T- 
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_ irtiSHAOf.'o.oi  CO  TO  tit 

call  s-:iiTTst,9-  st.c«  s:.7:«';t.,>c«st.0e'.st.«cisc.rc:sj iIClfmiRZ  t 

lK.PM2T<,0?tCl.LS.H, H.3.ISCS- ) 

CALL  COISLtCO.XCXSt t?C«S' ,T;c.«CC.Sl£,PHt,THE,CCiFF.GC> 

_CALL  P?(YG  aP) 

call  CLITGHCL.fCN.aMe.THt.SLfcCO'.TC  PRl  ' 

RM0FI‘*9M0  _ _ 

call  Sl(0?t$liAl  Sl'.A^  SI, SI  >1, Cl  SI.OXISI  .BilSlVLlISliLZlSltRilS" 
_ ll,9?ISl.SNISl.PH?T<.e?  Y<,Oj  S'.,F£  St.HT.Pl  Tx  ,Ph1  TX,  ICLFH , LI  TX.L 
XZ  TX,tTTSE,'l,A,;.t>,1.xTlSl.»CISl.LS) 

call  SP(VT<S’.VSHSZjArKS:?,OhZT<,  a?  lK.02ISl.8i  SE,A1  Sl.AZ  SI. OX  S 

it. SI  si.  ITts£  .y'-.h;l.h»  ,tl  S 1 . 3* I S 1 , SilSl . xCISi  , 2C I S£ .LZ  T 

_ ZK.LlISl.L2ISl.PlISl.P2Iit  .a;LSC.r(aSl.»C*SE.D,LS.AM.NX,NM.HYl,LX,S 
SM.VCO.Pl  TK.AGOSZ. AyasZ. ATCi2.A:KCZ.AACS2,ACCS2,ISGSE,VTISl.ACISl. 

XACRSZ.ATISl.AT^SZ.PSfiSZ.xCHSZ.xTKSZ.ZCHSZ.ZTXSZ.L.LJ-JLX) 

OVP- (VFhSS-VC-<S?1  / < (ACH/PKl-apiTl 

__  PRAT  1 «0. 0 1 ( t,  0 , (PCH/oTKI  I I ' 

CALCULAft  OVTk 

call  HClMX.PRATt 

hy«mx»hti 

call  T<<IS-(TK,P^TK,RZ_fK,Rl  TK.PhITK.LI  T Kj^L  Z_T  K,_HY  , a i EjL  ,L_S_» 

ISHAPitlSMTK  " ' . ' 

ir (ISMAgf. -o.oi  CO  TO  189  

CALL  SEdTYSi.li  Si.OY  Si  . JCXS- , XCXSE.  DEL  St . XCISE  , 2C  ISi  , ICLFM.RZ  I 

tK,fH2TK,C2ISl.Li,M,S.r.!;OSil_  _ 

CALL  C0(SL-.:0,  XSXSi  .ZCXSi  ,Y;c.YC:G.Slf  .PHE.THE.CC.PF.CCr 

call  PA«YC  PR)  _ ^ 

call  CLIYORCL.l'CN.PHe.THL.rLACO.Vo  PR) 

RHOFP.AHO  

CALL  SlOZVsi.Al  Sl,A>  SI.  St  SI. Cl  Sl.OxlSl'.B- ISl.LlISl  .LZISl.RlIS 

ll,0ZISt,SH:St,»HZ;K,»2  T<,3<  SE.Bi  SE,HY,pi  TK.PH1TK.ICLFP.L1  TX.L 

*Z  TK.ItYSC.O.A.S.N.R.VTISl.ACISl.LS) 

call  SFtVT<3’,v:MS? .ACHS?,PH2t<,r2  TK. 02151, BX  SE.Al  SI.A2  SI. OX  S 

it. SI  SI.  ItYSX.YCHSL.HY  ,ci  Sl.OYltl.SEtSl.XC.ISi.ZCISE.LZ  T 

?K,LirSl.L2tSl.RlISl  tP?I3l  .0-LS£.X);!Sl.YC»SX.a.LS.AH.HX.NH.HYI.LX.5 

Jm.VCC.RI  r<,A<-.P32.ATA'S2,ATC32.AC«.S2.AACS2.ACl.'S2,I3&SE,VTISi.ACI5l, 

XACRSZ,ATIS!.JTRS2.PFAS?jXCHSZtXTKSZ.2CHS2.2TKS2jL..Ll.J.<) 

OVK* fVT<5Z-VT<S)/. 01 

VTK=VTkS2  

199  CONTINUE 

C CCNVERT  RAOIANS  TO  OEORCeS 

'thctae.thEiaE/raoian 

_PMIE’»hU/RAJIA-< 

OTMXT AsOTHTTA/RAOi AN 

OPHIX  rCPMlE/RAOtjl^ ^ 

SIF«SVt/RAlIAN  ’ 

RETURN  A;N', . _ 

' 'USRCUIJNC  sitlTYP.RXTA.aCLX. 2CX ,YCX .DELTA ,XCHI .zchi.Tcall.rz.phiz 

X.OZI.LS.M.R.O.ISXG!  

C blVISIO'i  0'  TH£  TRUNK  INTO  SXCNXNTS 

PEAL  LS  _ _ _ _ 

dimension  2:K(32)  .XCX(32)  .ITY?*3Z).OfLTA«32)  .ISE&«32).XCHH3Z)  ,2CH 

XKSZ) 

DATA  PI/3. i<lS926S3F 

NST0p=32  

C IF  FIPST  CALL.OONPUTE  PARTIAL  TERNS  AND  NunaEP  SECHEHTS 

IF(ICAlL)  ?0. 30,20 _ _ 

30  kLSM.n.?»LS 

C,  beta  is  ruRKEP  segment  arc  ancle  

eETArPI/2,/»-LOAT(UI 

C OELX  IS  STPAlCNf  SEGMENT  LENGTH  

OCLX.LS/FLliTtZ'N) 

8CTA?=).  3 3Hl»SINnPTA/2.)/3ETA 

C NOMRlF.ING  of  SX'.NiNlS  AOCOROIMG  TO  ThEIR  position  in  The  trunk 
00  11  I'l.NSTOP 
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k 


IFd.LC.NI  ['^^Gdl  il 

IFd.CT.N.  J'd.  IGE':d»*^  . _ _ 

IFd.CT.N*  1.1>l3.I.L".H*?»>*dS'Cd»»l  • • 

iFd.GT.N*  i.iF.J'tN.HjdfEGdm. 

lFd.GT,^.*(N*•1>.A^I^.I.I.t.  IS£Gd)>S  

lFd.C1.3*'(*3*'i.*x"'.:.L:.?*«lcinI3£Gdl=G 

lFd.GT.3*lN»‘il.iHa.I,t£.3»**«<i*1l  IS‘G.<  I >«r  ^ 

lFd.Gt.3*'l»i.»M.»Ma.I.Li.<.*IN»HII  ISECd>:( 

11  CONTINUE 

c evaluating  »<»o>:^Tt£s  o'  segnemts 

C 1TY»*1  FCR  CURVES  S;GMENf..O  FOR  STRAIGHT  ScCHENT 
,e  *CX  AND  ?C*  AR-JI  AV>J_E  COORDINATES  FES®.  OF  The  SEGMENT  CENTER 
C X?HI  ANO  ECHI  IRE  t ANO  ? COORSiNAtiS  »-;SF.  of'  THE  CUSMfON 
C PRESSURE  CENTER  FOR  A segment. AHEM  I’  IS  OUT  0^  CROUMD  CONTACT 
C DELTA  IS  S-CMENT  CENTER  ANCLE  RELATIVE  TO  CG  ' ' 

.EO  CONTINUE 

oefo.sfo^fifsiniphiei  - - 

_ 00  10  l«l,NSTOR 

RCO«lS£Gdl'  ‘ 

CO  TO  <1.2.3.A.5.6.7.g>,  KGO 
C CURVED  segment  “ 

C IF  NOT  INITIAL  CALL  SXIR  CALCULATIONS 
1 IFdCALL)  R.IOO.R  ■■  ■ 

_100  ITT.-!H«1  ~ 

OELTAdl*rFLOATd-ll*0.5l*3Et4  ' 

COSOEL»COS(3-LTAdl  I 

XCXdlF- <R.jHH0>I‘COSQ£L» 

ZCXdtcOEIdlNCfLTAOT 

XCHI  (Ii«.  (RLSH*0EI  •8£TAE»C'OSOEL)' 

2CMld|.IC«dMlETA2 

60  TO  9"  “ 

c straight  SECHENT 

i ~ ITYPdl.C  

XCxd)«-RL3H.<FL0AT  (I>1<NI*0.S>*CELX 

ZCXdHOE 

xCHidiixcTdi  

zcHi  di  *ic<(  iVroTs  — — — 

CO  TO  9 

C straight  SEGHEnT  ■ 


_3 ITTPdl.O  

XCX  (I  I •(FLOAT  < I-N-N-l)*0.5'l»OELX 

ECXddOZ 

XCHidIrxCXCI 

zcHidi.Ecx<ji«o.F  

CO  TO  9 

C CURVLC  segment  . 

C IF  NOT  INITIAL  CALL  SKIP  CALCULATIONS 


,>  . IFdCALLt  1..00.9_  _ _ 

<•01  OSLT  A ( II  »,(F1.0AT  ( I-N-E«*1-1  t ♦O.S'l'PErA 

SlNOEL«SIN(OELTA(II  I 

XCXdIsRLSR»3ZI»SIN0EL  ' 

zcxddozi*:os(3:LTAdt»  

XChI  dl  ■PL;M»0ZI*nETAZ*SlH3EL 

ZCHldMZC<dJ.*SElAZ 

60  TO  9 

C CURVtP  SEGMENT  

c IF  not  initial  call  skip  calculation: 

* IFdCALLI  1,500,9 

500  ITvPdi.i 

delta  (I I . (FLOATd-Z  •N-;»M-t>«0.5>*BETA 

tOSOFL,COS(3"LTA(Ii  I 
XCx  d I =RLSh,3tI*COSO-:l 
ZCX  (I  li-o?  fSlNIOEL  T A(I  I I 

xCmI  ddRlSH.nzI'COSOrL'RFTAZ 
ZCHldltlCdll*  lEIAZ 
CO  TO  9 
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C STI»»ICmT  SEOtlENT 

» ITTO«I»»0  . _ _ 

iCXCI»»-32  j 

*CHm)*xc<tr» 

?CMt(ii«?c»tii»a.5 

' GO  TO  9 

C STR»ICHT  S.ECf'tHT  _ _ 

7 iTYp»ri«o 

xcx  tII*-»fL04TTI-3*N-3*K-l»*0.5*»0.'L*_ 

7CX(Il«-07 

xcHHH  = xc<«t) ; 

7CHI t II»2C<<: I *0.  f 

GO  TO  9 

C CORViO  SiGHiNT 

C If  HOT  INITIAL  CALL  SXl®  CALCULATtt»l5_ 

> IfCCALLI  9, too. 9 

ITVPdI.l 

DELTA  (tl r (TlOA  T ( T - 3 •N-L •«- 1 1 4 0. S » • BE TA 

£IN?EL«SlN(0*LTA(n  > 

XCX (1)«- (P.SH^OII'SINOELI 

2CX(Il=-0ZI4C0S(aiLTA«T» 1 _ ^ 

XCHI  (I»«-  (TLSH*0214SIN0- L*9rTAZI 

zcHi  (i>  = 2c<(i  I ‘^ctaz ; 

9 CONTINUE 

l6_  continue  - 

RfTURN  ~ 

tNO  _ _ _ 

S099C0IINE  ’kdSHAOt  HIl.Ll.  LZiMY.A.e.L.LS) 

C TRUNK  geometry  OALCU'-ATIONS 

' real  i.tL‘r,LZVL"S  ■ 

aTOL*.l_  ^ ; 

l®«HY.Li. 0. Jl  GO  TO  111 


tIERATICK  POR  TZ 

_COMOU1E  INNiR  RADIUS  Of  CURVATURE 

SZfSORT  IA*A*0.  Z;  »HY4HYi' 

YtERATION  "loop  fo"fi'LZ.GU,oT7RZ 

CO  lOZ  1*1.50 

PmIZ«A85(  A IDS  ( A>*4X1<-1.  O.AMIKt  it.  0,  ( IRZ-MYl/RZ)  M > I 

_ SIf|PMZ*SIH<OMlZI 

COMPUTE  CUTES  RADIOS  0®  CURVATURE 

.M«(l  A-RZ*!! -iPM?!**  ?■.«'».  MY  4(04  MY  

PHI  1«APS(  AIDS  ( AMAXl  (-1.  0 , AilNKl.  0,  i 4R1- HY>  8 ) /R  T)  till 

xs^A-RZ’siHPHZ  ; 

IF  (XS.LS.O.  01  PMlir6.Z«3n'5Z-®HIl 

LZ»L-PMI  1 ‘RL  

RZS  IE  resultant  radius  for  computed  LZlIN*ITERATlbN 

If  (ARSCPHITI.LT.  l.OE-Z»  PHtJ._t.Of-Z 

PZS*LZ/PHIZ 

TEST  If  toleranc:  .gt,  erpor : 

If (AMS IRZ-RZS t .LE. PTPLI  GD  TO  >0 

RZ=  (5Z«RZS)iO_i5 

lOZ  CONTINUE 


iterated  50  TIIES  MITmOuT  SUCOESS. ERROR  RETURN 

_111  CONTINUE  _ 

NRITE (t.9011t 

_9011  fO»MAT(nx,^IHffASAnLE  T^OHX  ftOMETRY  _?//» 
ISMAPEiO 

friUAN  _ 

C trunk  OK.RETURH 

50  ll*L-LZ 

ItMAPE.l 

RETURN 

n 
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ENO 

SU0»00TINt  .DZ32‘<i'I»D*AH»Br.BETAD?I.llI.L;i.»lI,R2I 

X.SlNPMRt.PHIJ.RZ.rillJ.BEtA.KV.Pl.PMtl.ICALLtLlfLS.ITVP.OiA.E.N.N.V 

XlKl.ACHl  .lil  

C iHITlAL  ASSiSS'(iST  OF  AOiAS.tfPLllMtS  ASSUMING 

C NO  GRCUNO  COnTJ'T  , _ 

RIAI.  Ll.L2.l.lItt.2t.LS 

OtliNSlDN  tTT»«!?l,AT<I(32t.VT«I(J2l,ACHlC32l  

NST0P*32  ’ JB=E 

C COMPUTE  GfOMETTY  TERMS  

SINPH2«SIN(0HI2»  ■'  " 

SINPHP«SlNaH2*5^ 2 

02=0/2. tSINPKR 

U02«0YL«*02  

BC02*BErA*32*02*0. S 

*•3* IA-SIR®M=I /(9 

C COMPUTE  AREAS  If  TRUN)  SECTORS 

A1.PMI2/2. 1*R2»*2 ~ 

' A2=«»2-HT|/2.0»SIUPHR 

_ A3«PHll/2. 3*R1*‘.2 

AA«X*B/2.a 

AE«  «A-StuF-IR-*l /2.  0 • _ ^ 

XI.SInPmR...  0*  (SIU<phT2/2.  on  ••2*R2/<3.0*PHI2» 

J(2«a.  S6S67*S.I*iphp 

X3*SINphp.4.,o»  is  In  (PhU/z.oi  »••^•Pl/l3,0♦PHIl> 

*A«A.a.!33I35*<  ^ ___ 

YS*SIN»mR«Q,  3 333  33*(A.*S:»lPHP.-i» 

S«2.0»L>»e. 21313*02 

AAi;Al.A3«Ai-A2*A<. 

AT«A1  ♦X1.A2»<2»A3*Y3-A<»XA»A5*«S ^ 

IFIlCAtL.  CT.  01  SO  T0"23 

C SAVE  TRUNK  OEOIETRYjrRMS  FOR  EnO  TRUNK  CAtCUCATtONS 

«ll»»l 

_R2I'*R2 : 

PHllItPHll 

_PHr2l.^HI2 

L1I=L1 

n2I«E2  

Atl^Al 

_ A2I=A2  

SlNP“2I  = SI'|eH2 

_SIHPhRI.SIyP^R 

X1I=»1 

X2I.1IZ  ^ 

AlMA2I=Al>A2 

_ 021-02 ^ 

$I«S 

8ETACI2I»32IA*n2  

X12=(Y1*A1.y2*A2»>AIMA2I  " 

_ DRAMA  R1tIO*0,E«<121  ♦A]HA2I»RETA 

t202>‘ni  = 02*02*0.  5‘OETA 

_ 002I*0£LX*32  

e002I-:e£TA*02*02*0.5 

_20  CONTINU-  _ 

C COMPUTE  TRUNK  '.iSME-lT  A RE  A , VOIUME  . CUSHION  AREA 

00  103  I*l,NSro“  I 

ifiitypiii  .eo.  nco  to  lii 

C straight  part  if  trunk  _ 

131  ATKIlInAA 

VTKtIIMOrtX*AIKIHI ^ 

ACMIlllOO? 

CO  TO  102  _ _ 

C CURVLO  PART  OF  TRUNK 

112  1F( ICAIL.GT.OI  GO  TO  102  _ 

ATKIIIKAA 

XE'AX/ATKnil 
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WTiCKIt.OE  rA*lO/?.»»E  l*»TKl(I> 

»Ch1(I(=b03E  _ 

102  CONTINUE 

103  COMTInU-  _ _ 

KfToPH  S£N3 

SOOPOUIINE  hCEE.KI  _ 

c soeiooTiN-  10  lALCULATt  positions  oe 

C SJO;  T»UN<  Loa-S  . 

■c  The  position  «-:tPR'isEO  pt  nyi  cepenos  on  pRESsuats 

C I.E.  H»/MVI  = r t^OH/PKi  

C FORCE  INPUT  PRESSURE  RATIO  'lE’NsEM  0.0  4NO  1.0 

AH0*1. 001  tlHl.-.aSS  }AHl;-.l,03 

XpAHIni 11. 0. Anat 1 I 0. 0 . >1  I 

_ 2.AH0*AM1 ‘X* AH2* «• X*AM3*X*<^X  

IF(7.Lt.O.  l»E  = 0.l 

_ IF  12. ot. 1.01  zn.a 

return 

eno  _ 

SU9R0UTINC  CL  <YSH,  KLN.OHtE.TPiTAE.SLA.YCl 

OlMiNSIQN  ;L-.t32  I.YCIJ’  I.YOhTJ?  ) 

NSToF=32 

.c  calculation  of  trunx  groumo. cl-aranse. for  each  sechent 

ICLNStlCLN 

cosrcssCosiPNisi •cosi.lHilA.ii 

C CALCULATE  SE0N-:nt'  C4’ 

00  161  I.l.NSTOP  . 

YCMIII  jSLL(It*»i(II  ‘COSCOS’ 

c _1F  NEGATIVE  SET  GAP  TO  TRRO : 

YGHdlsANAXKYOHd)  .0.01 

IFIYCMdI.LE.O.OI  1CL*IXICLN«1 

IGl  CONTINUE 

._XCLN»ICLNS : 

RETURN 

ENO  _ - 

SUBROUTINE  C0(SLA.«CX  ,Y:x,YI6.XCG.S:E.PMl£,lH£TAE.CC,Fr,CC) 

c THIS  suosoutin:  calculatts  x ano  2 ccoroinates  of_tm£_grouno 

'c^’roimt  corresponding  to  each  SEGnEnT.  for  a paPIICOLAP.  acls 

C_ORI£nTATION  

OIN-NSION 'SLLfs?  'l . YCXd2  I.ZCXdE  l,XG<32  1.2CI32  » 

NST0P  = I2  _ 

C'CALL  EHATRIX  for  SpaCIAL  TRANSFORMATION 

C 9HATRIX  TRANSFORMS  A VECTOR  FPQN  VEHICLE  FREAHE  TG  INERTIAL  FRAME 

"C  - - 

_C_  _ . . 

C CALCULATE  TSAl.SCENOENTALS 

csiE.ccs(si:i  

CPMlEsCOS I’HIEI 

CThETAEsCCSIThETAEJ ^ 

' SSIEiSINlSIE) 

SPHIE=SIN (PMIEI  

STHCTA£»SIR(lHEIA'el 

C . 

C COMPUTE  translation  MATRIX  ELEMENTS 

eH=CSIE*C’HlC  ♦STMiTAETSPHlElSSIf 

ei?tsPMi£*:THEiA-: 

_ ei3s-SEIF*;aHTE.sTHET£:*SFM:ETCSIE  

821»-SPHIE*CSrr«SSI' 'CPHITTSTIETAE 

BEExCRHlCTCTHilAE  

B?3iSOHU*  SSIEtCSIE  •rFHI£*srM£TAt 

e3t=ssi-;TCTM-;TA: ^ 

b32«-STHi;TA£ 

a33  = r,siETcrH£TA-:  

C 

c 00  LOOP  OF  All  segments  to  grouno  position  _ 

DO  tOL  I.l.NSTOP 

*cxcc=ixcx(ll-cct 

T' 
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ZCxrf«IZC»III-F^I 

C _ . 

C C»tCUl»Ti  VtCTD»  0«  F0»  $rG1€N» 

simi  t(YC‘.««f.xcc*ii2»7c»rr*6<2» /02?-Gc  _ _ 

SU‘.GG*<Sl.i.(It»CGi 

_C  ■ 

C CALCOLOi  K-G«1UH3  C00“.0lNATi 

»c(:)«^c*c';»»ii-sL‘*GC*jji»r';xxr*p3i«xcc  

c 

C C»LCUL*r£  3-CB3'J-iO  GOOOlIMiTE  

2C(1»  *«C*CG»  ll3-S;.<.G0*321*2C*ff*833 

ItU  CONTINUE £ 

■return 

two 

SUSROoriNE  PRCTOl 

c USER  sPECiri£o  CROon  rroxile.  _ _ 

'c  CLEVtTICN  YG<r)  :s  -yoRESSCO  as  t rutiCTlON  Of  X RNO  2 COOftOiNXTRES 

C ox  GROUN3  POIKY  I.I.E.  »G(I1  t <0  TGI  II  ; 

OlNiNSIQN  YG(32  t 

NGTOe«S2  

00  105  l*l.NSr0P 

_C  SET  rC*.  RL»T  TtyAIN ^ 

YC(  Il’O. 

14? CONY  l_Nui__ ; 

RETURN 

iNO  , 

SuaRCUIIN?  ■s'»<vT<,i/:H,«C«,om>,^2,22I,8£TA  ,A1,*2,0ELX,SI,  ITTPTt?N, 
_ XHr,O»03-<5:.3<»N»q;  ,«t'A32T,XC'*I,7C“!.l?.LlI.U2I.Pi:.R2I.0?LT4,SINP_ 

' XHei.XCXtO.LG.AM.NX.NK.HYI.'.XtSN,  VCMO.ftl.A0AP,iT<AT,*rKCH,*Ti<CN,ATK 

XATC.AUCMC.lSrG.XT<I_,  ACnI,a:h»,  AJ  xC  N.I  . A T K C NR  J f £ R.I . XCH,  X T<.  TCH , ETKt 

XL. Ill'  

C CALCULATION  Of  AREAS  _AN!)  VCm4i5_ASSOClATEO  nITm  ACLS.XNOWIHG  ITS 

e ORIENTATION 

REAL  LT.LS.LY.LlI.LTl.L.Ll  _ 

"OINENSION  YT<09(32  t.VTKRAljZ  ),VT<R5(32  I.PERIUZ  ),XCNt32  I.ZCNI 

■ X32  J.AICCNldT  (,AT<C'iR<32J.4TRATH32  I.AIXATRUf  I .VChI  U2  >,A6A 

XP1(32  I,AG1PR«32  r.'xYKI32  ».2Ti(<32  I..MXP.(32  1 . ACHI  ( 32>  > XCmI  « 32  )i 
X2CHI(I2  l.Y5H<32  l.TffOI  J2I.VTXH  32t.XCYt32).ITYP(32).D-LTA(32l. 

■ XACNRt  321  , v:HR«32t  . Y T K(  r 32  > . A KCMR  ( 32  » . ATkCmI  ( 32) 

_ _ NSI0»*32 SPI»3.  1I.1S9ZG53  

AT0T4L».6 

_C  CONPUTE  partial  TERNS 

SINPnZsSIN (“HIZ) 

SIHPmR  = SInRn?»R2 ^ 

02»0’  0.5*S  INPNR 

l.2IS0«021*32I 

A1NA2.AI-42 

6tTA2.1.  3 33  3 I»SIN(  RET  A22.  Y/aETA E 

RISH.LS'O. ;o 

_ aos»ah»oelx/si  

Ae3SI>AH.pETA»32I/SI 

fcHN«»LC4T 



C PART  I I value  ox  v:h  ano  agap  _ _ 

'c*»»* 

00  IT  I.I.NSTOP  - 

C TEST  TOR  TRUNK  SEGMENT.  WHETHER  CUNVtO  OR  STRAIGHT 

iniTYPtli.EQ.tl&O  TO  ll  _ _ _ 

C 

C STRAIGHT  PART  _TX_T PUR*! 

"c  ■ ■ 

C calculate  Cushion  SEGMENT  INITIAL  VOLUME  _ _ 

13  VCh:  ( I).  ( Y -.•u  : I "02- AINA2I  'O-LX 

C CALCULATE  SEGMENT  GYP  AREA  _ 

AGAP]  1 1 Y « ( YCH< I ) -HYI "OELA 

GO  TO  10  • 
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c 

C CUtVCO  or  T<UN<  ^ . 

c 

11  wchi (It •V6MI I 

4CAPI(1I< IVOH(I)«HVtt«1£T4};i 

* ••••••*••*••" 

c PART  i R VALUC  CALC'JIAT  IONS 



•e 

e trsT  roR  croon:  contact  at  iacn  sicnint  _ 

It  contihui 

^_rORCl  VOLONS  ARCAS  .Ct.O 

vCHi  (Iiianki  1 0.  b iv'chItii  > 

ACARI  l|l•ASA<(l  (0.4.  ACAPtItll  _ _ 

‘c  TEST  SEC“iNT  rcR  CONTACT 

_ JPdTVPdl.lO.  l.AlO.tCMdl.lE.MTll  CO  TO  1% 

■ irdTTPdl.SO.O.AMO.TCMdl.LE.HTl  GOTO  IJ 

c _ ... 

C NO  CROUNO  CONTACT 

c SET  CONTACT  AN:_RJ*1OV'_TER*tS_T0  iERO 

ATKRdOO.O 

ACHRdXO.] 

VTNRdXO.O 

VCHRdHO.J 

ACAPN(i>*0.3 

ATKCNKK-O.O 

atkcnriii -no  ■ 

ATKCNR(tt>4.a ; 1 

~ ATKATR(I»«1.0 

PERtdl'O.I  _ _________  

c set 'oisTANCts  tTr  fo  r^ti  f*oM«  viruTi 

»CHd)«xcHtdi ■■ 

ECNdIdCHKIt 

2TKd>«:cH(tt 

XTKdlOCHdl 

»0 T E T»ONK-:OSMIOM  - ATMQSPwftt  8Att  Q API AJ| 

iPdiTOdii  rt.itTit” 

_1*  CONTINUE  . 

C NO  CONTACT  STRKCNT  SICfloRl 

ATi(CNl(ll»'LOAT(iriX((UI-tX»/JM«1.0»^HI»APS  

ATRAT  I(II»TNN‘A:S-ATXCHtd» 

60  TC  IT  

It  CONTINUl' 

C HO  CONTACT  CURVlO  SECTIONS  

ATRC-KII  •‘AOAT  (iriX(  (Li:-LXI/SH»l,  0»*NH)»A§6SI 

ATNATI(n.TNII*A10Ui*T5CMIJU 

CO  TO  IT 

• »»»»* 

c-'  ■ ' ' 

C TRUN<  CROONO  CONTACJ 

e' 

C CORViO  PART  or  TRIJNX  _ _ 

C CAICULATC  O'roiNATlON  ANCl'S  rOR  SEONCNT 

1% PHI  j»Acos  ( (r:i  - (MTi-»cN(n»  t/Rio 

PNlAfACOS  1 lRlI>(Nrl>TCNd  l»  I/Rlt  I 

SlNPMJ.SlNIRHi:!  

SINRMt,>SlN(PHt  41  - - 

C COM»01f  PAATIAt  TIRNS  ,■  

OPSP*  (OII-RTI*SlNrHII  

(-R5PT"0R5P*:RSP 

COtO{(.*C-!(:ELTl(il » 

_ siNect»siN(jiLTAdn 

Of  n»i(ii4s  ia:*otsp*5Inoii. 

UrORCS>RCT  IT'ORSR'COSOEA.  

C C0NPU1E  PjPOVAL  C':CT0'’3 
At>Rf|«NII*PHt1>0.t 
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I 

I 

I 


t 


{ 

! 


f 


f 




*<•  I •0.5*HlI*$lM«»Hfc  _ 

A10»i« 

-- 

C COH'Utt  I!CT0»  C»»lT»OIO$ 

m^SlNPMPT-l.  J313S  J*  (‘IMlPHl  J*0.  »*••?)  ••21/PmIJ 

*r»iiNPMP:-o.  3un3*P2:*siHPMi 

»»«SINPMP1»1.  3333  331*(SIn(»KI<.»0.  5»  ••21  •P1I/PMII| 

K*"1INPH«1  •0.  3133  33  I'PlI'SINPMl. 



xii*«a 

pii2«pr»fl.3 

IP«PHI<..13.PI12»  CO  10  30 

c tP  PMifc  cp!*t23_than  »o  otcpEc^j.itLJlo.j.o  oecpcci 

Pm1i..»II2 

tTNPHk>3IN(PHIkl ^ 

’ *t0ttPlI-HT:*T6M«nt»PlI 

*10"1INPN»I.0.3*P1I 

...  . xii'SiKPHP.i,  333111*  (?ifi(.p>i:kj*..#.ij;2»r*n/ptiis 

• 0 ccMTiNur 

JL  tWi»uT;  .T»LN<  ^ 

iTKP(tMA614'.A3.A« 



■ *rP»<»o‘«e-***<’.AO*3(-*3»«3i/*T((pa» 

*CP3««****»-*2*«?‘«lJ*’><»l*HO*»10»/»TKP«IIj 

'c  coMPuTi  r3uM<  vouun:  cHtNce 

Vtl<Pt(tl.<!»*4»(O.0.3«»;PI»4f<t<H 

VTKP*  {IU4f«P(II*3tt4*(d4n.3*MP«l 

pTXPdHl.'VTPPKtCVTKPOXIJ 

l(CP4(4t4l»-4r*(7»/46M42 

C ^OUPUTC  TPUN3  {3IT  40!43  , 

iTKOHl  (II  ••4,043  (I»  IK  ( (L21-U»-«2I*PMI3*/SM.I,  01  •HH>  •48031 
• 4TK4tT(tl  ••1043  (I  P I K ( (Ul  l.LjlK.»t,04  T (WK»  1 1 »lwml»PMlH /fH  . t,OI  »NM 

11*49031 

ITKCKPI 1 1 ••1.043(Ip:3«  (121-13  1 /3*<*1,ci*NH1»4005I»4IKCMI  (^I 

*3I(I3P(II  •3»C(*4333I-4tKCHI«II-43K43I  (II-43kCM»(3I 

PrPl<tl«9Cr4««OK3».02!.’(lI»31’«PH4l 

C‘e6"»U3£  CON34cr  »tPI3!3?P 

43KCNI  (XI  •1£T4»1.«*  (02I3Q-OP3»2> ^ 

'c  eOH'W3f  3KU‘(K  :0'l34Cf  4«4 

„ AIKCNP  ( 1 1 • 3’.  T 4 • 0.  9T.(  (n2l*»lIiJJt!tP!ikl”.2rP,WiQJ 

*CHPI|I>43<:nI(1I 

0 COf1»U3£  CU1HI03  VOLUXt  CMAHCt 

»CM«(II*-<!i34*41147»  (0‘O.J.KCPI 

C COKPU3C  C4»  4P!4  CHJsy; 

4r,4Pf  (jn4;iPi(ii 

c oii3»((Ct  ()*■  si  iUNr  ppisiuPt  CfNliM.  *^0”  cosHipH  ccNTej 

2f  lip.oil.pil  MIHPHk 

AP1"C21-P2I*3IM«>m.»  _ 

*»2»1.3Jlll3*S|*(((llTl»0.JI/lC34»(PP**i-(lPl*»j|/(PP*«P-P«l’mil 


»C0*I3I0(II  

"■  CO  30  (01 ,11,2l<Ak,C5.;i,21.C4ltKCO 

Cl  *';m(ii..»i,i-(-'iiopci_ 

2CM(II>9r033N 

»3K(I|..«l.lH.KK2»C030[t  

23K(ZI>K(2*3lNO(l 

CO  30  IT  

CC  KCHII (•PlSH.If 0«|N 

7CH(l|.'1IDIC3  _ 

»3K(1  l•^lI■I.KK^•1IN0^C 

23<(Il.««2*C030{t  _ . i 

CO  30  IT  I 

Cf  KCHlIKKlCXfieoKCS  ] 
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i 


iCHdifatonH 

2rK(iifxx;*stH:;L 

GO  TO  17 

at  XCM(Ilfi«LlH‘>>IC3?SM 

XT<(  IM-iSlJ-^-xxJ^StNOE'L 
2T<(Il.-xx':*CCS3rL 
CO  to  17 


tauNX  CSCUNO  CIHMCr g 

StRAiCHt  PXRI  jr  TRURX 

tl  CONtlNU- 

C COHPbIC  OiroRMinoN  XNGLCS 

RHt«  ( (97- dY-VGHd)  ) l/R2t  

RMIlclCOS (AIXYl (-1. O.AnlHl (t-O.RHVI 1 1 

_ RMY*  ( (9t- (HY-TOMd  ) 1 1/Rll ■ 

' f>MIi.aXCOS  I tltXl  (-1. 0,  CHIMKl.Q.RHYri  ) 

C 00  TRXNSC-ROENtiLS  0»Lt_0MCS 

SINPHl»SIN(OHl Jl 

SINOh..sIn(?hI(.I 

C COROUIE  pxSTIAl  TiR-S 

DPSPa(0?-«1»Sl>(»M3l  

C‘>SP7’»09SP*0’SP 

coscit*coso-:i.T»d)t : 

SINaiU'SlNOiLTAdl  t 

e!09SNTorrA’*0RSP*?INj£L 

BSC9CS»9iT»7*09S»*COS0eL 

C C0f1«UTi  Pi«OVt'.  SJCTORS 

At»R?*97»<>HIJ*a.5 

*7«(R7-hy  «ygh(I) I»3««»R7»5INPw3 

X6»(A7ta6*X7 

_ 

■ A9«(»l-HT  *YGHdl)»0.5»Rl*SINPMG 

_A10»A9 

AU«A« 

C COMPUTE  S'CTOP  CiMTROlOS  

xeaSIuPMA-l,  353dl»  (SIMIRmFi*  0.5  » ••7>»R77PMI3 

X7*SINPHS-).  J33333‘P,7*SrM®M3  _ _ ' 

X7a?INBMP«l.  33  33333  •(SI;i(P'(Ii.*0,f»»*7»*Ri7PMlA 

X9«5INPHR.J.  333333  3*RfSINPHfc ; 

“ XlfiaXO 

xn.xj  ^ 

PII7«»I»0.f 

IP(PmI»,lT,piI7I  CO  I0_70  

C If  pmia  is  GRClTiR  IH4N  ao  OiGRiES,  SET  TO  90  DECREES 

__  PHItaPIIZ  

SINPM,.SIN(“HIV1 

_ AlOaCd-MY.YOMdl  IBSl 

X10»SIN»Me»0.5»Al 

Atl»»l»»l»»MIl.»0.  • _ 

X11«SINPHP.1.3  333  »3  3»(SiM’«P>1IA»0.5»’»73*Rl/PHIA 

70  CONTINUE  ^ ^ 

C COMPOTE  TRUNX  IREA  IHANC'" 

ATXPdl  .Af  NATaAl-Aa 

AIXPEdHAiTATYAl  I-AIO 

C COMPUTE  taiNX  V3LUM-  CHANGE  

VIXAAdI  tATXR(  II'OCU 

VtXB'Tdl  lATSRTdl  ■'Tttr 

VTAKII  .7.  •VTXRA  d I -VTXRO'Vll 

C COMPUTE  CUSHIOI  VOIJHE  Change  _______ 

VCHPdt.-Oax’ASHAT 
e COMPUTE  TRUNX  -YIT  AREAS 

ATXATIdl  iBLOAT  (irlY  («IJ.L»(.X»FLCAT1NX-I)»SM-R1»PMII.>/SH.  l.Ot*NM» 

1*A0S 


I 


I 


1 


! -1 


1 


i 

1 

i 

I 

i 
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*TKCMHl»*'LOtl(lf  I»l  a?-L*-R?*»Ml31/SM»i.  )»hM»«AOS 

*'(.o»T  ti^i*  I a z-t* !/■;«.  i.oi  •NH)«Aos-ATKCHi(t) 
*T<»l«tIl»?NN*«T>-*TKCHr(I».tTK»TI(ri-»IltCM»II» 
pt«ni»*z.  *3^1* 

C'OHPUT;  I’UMK  :3MTiCt  ‘•'SlrtiTf* 

COHPU’E  trunk  CINTACT  *»£» 


ATKCnIII»»T?«SI*(OhS*C'L« 

*TKCN®lll*Tl*5lNPH4*CEt» [ 

C COHPO-ie  GA®  ASIA  chance 

ACHRtI)»AT<CNIIII  

. ACAPRIIloAlAPIltl 

KG0>ISEC(1I  ^ 

C CONPUTE  SrCNEMt  CONTACT  CENTER  Qr  PRESSURE  fOR  CUSHION  AND  TRUNK 

CO  TO  <17. >2.62.1T tlT,A6,A6.l7>|KCO  

62  XCH(I»«(C>(tl 

_ 2CH(n.a.J*02-R2‘SINPHS>  

XTK(I>>XCX(II 

. 27k  I It  «0  24  l.S»<Rt»STNPH<.-o.2«SIHPHH ; 

CO  TO  17 

66  XCKIU.XCXIII  _ 

2CM(lt«-0.}»<02  R24SlN»H3i 

XTKTl  IttcxdT  _ _ . ' 

2TK<1  l»-  ;0240.S4(R143IN»h«.-A24SINPhI>» 

17  CONTINUE 


PART  3 SUHMArnN  OP  ScCHTHT  areas  VDtUHE* 


SET  total  area  ASC  TOkUNES  TO  ZERO 

ATKCN»0.0 

V’K»0.0 
ACHiO.O. 


ATKCH«0.0 
>TKlT«0,g 
VCHiO.O 
ACAPiO.O 


ATKATC>0« 

ATKCHC»0.  

C LOOP  CN  3EGNENTS  ''0  RlNO  TOTALS  OP  AREAS  AND  VOLUMES. 
_ 00  JO  l.l.NSTOB  ^ 


VTK.VTK4  (VTKKI)  -VfKPlIlT 
.ACM. ACM.  (ACh:  <I) -ACHR(nj_ 
atkcm.atkch.atkcmh  i i 
ATKAT.ATKAT.ATKATl (II  _ 


30 


VCM.VCM4  (/CHKII-VCHPdTI 
ATKCNiATKCN.ATkCnI (I) .ATKCNRItl 

agap.acap.  (acapi <:i-acapriiii  ~ 
_atkitc4Atk»tc4aikatr(:) 

ATKCMC.ATKCHC.ATKCHRdl 

CONTINUE  

ACAP.AMAXl (ACAP.T.OI 

VTK.AMAdO.O.VrO  

VCM.AMA<1(3.0.(VCM4VCMOIT 
.VCM.AMAXI  d.O.VCHI 


ATKChsAMAx  1 10.  00  0. ATkCHI 

ATKAT.ANAXKO.  OOO.ATKATI 

aCH.AMAIld.O.ACHI 

ATKATC'AhaaHO.O.ATKATCI  _ 

ATKCmC.AHATKC.J.  ATkChCI 

.C_.  POACt  SUN  OP  N07ZL-:  AREAS  TO  Re  EOUAL  TO  TOTAL  NOZZLE  AREA 
SUH.ATKAT. IPKCM.A TKATC.ATKCHC 

ATKCH.AIKCH/SUH4AT0TAL  7 A T < AT . A TK A T/SUH* ATOT AL  

atkcmC’Aichc/sun'atotal  eatkato.atkaic/suh’atctal 


r 
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XYCC,CVCH«>,«C‘1.v;H.SINK»TtPHlE,0‘»Hl  ,THtT*E,OTHitA,CKi(.**T,&C,Mrl,C 

<tCf  nAE^.RPi,  v[i\,  sic,rho.a(;ap.ach,atk 

XAT.ATXCH,  ATXJN.ATXArc.AKCHC.CAC  .C£‘«rx,CCAP,  ACmR,4T<CHI  .AfXCN 

CA.P£RI.*CM,7CH,XTi<.?TiC,ACMl.t0h,C''LTK,QTXCH,0T<AT,0CHtTI  _ 

C OYHAHIC  »4K  version  EAR  EM4. 

C STAtt  EQUATIONS  POR  TMi  OV'lillC  system 

REAL  MASS  - - 

C POltCWtNC  SUSRTuriNIS  ARc  CALLE}  TO  UPDATE  VALUES  OF 
C EORCtSi TCRCU; S AMO  «LOWSt  GIVEN  Twi  NEM  VALUES  OF  THE 
C state  VARIAILC* 

OlMENSION  3£TT  (ITI  , ACHI  4)?t  .ACHTnZ)  .AKCKI  IS7)  .ATKCNROZI  .P£RI(S2~ 

. XI,>CP«32»  .ICMJTTI  .iaMJJj^jJ!Ti<(32»^,TGH4  J2» ' 

CC»- 1.17  jAATPN.;  796  " SCOT*. 6 TCTA»i<.  SQVENTiO. 

^CTC..A  JFE.O.  _ _ 

HOC'l.  IPHA<15S.  SUaO.OS 

_ A»LT<s.698  SDAMPC«3.2  SAIEAN«.eT2 

C£NE7«0.  SNST0P«32 

C_SVBR0UTINC  to  pimp  flow  ANQ  PP-:SSUR~  VALUES  DURING  OTWAMtC  SIMULATICN 
TIRMOO.  O /THO 

C PLEMUN  TO  tTUN<  FLOM  

SIGM^liO 

IE« (POLM-P’RI  .LT.  0.  0 t SICN»-1.0  . 

■ OPLTV*S1CN*CPT*LOLT<*SORT  «ATSITIFHO»<PPLM-PfRi I » 

C trunk  to  CUSM13M  flow  

* SICnVi.  ' 

1F(  (PTk-oc-O.LT.O.  ) SICN»-t.  

OT<CP«5I&N*CTC»30I1T  IA9S<TIRH0*«Pri<-PCHI ) ) • I AT<CH«  S.  1 6647*  ATXCmSI 
C trunk  to  ATH0S*HEP£  plov 

’ ■■  ■ SIGN»L,0 

IE(PTK,LT.9.0IStCN*-l.  ■ 

QTK4TsSlCN‘CrA*S1Tt  (TIRHO*A3S(PrK»J  »(AT'irA'T«07V46ir7'*ATKrt7rj 

C_COShion  to  atmpspher?  flow 

SICN*1.0  ■ 

IFIPCH.lT,  ).  01  S:GN*-t,9  

OCH4T«AiAP*CGAo»SORT  ITIPKO^ABS'iPCmI  » •SIGH 

^FORCES  ANO_TOF3UES _ASSOCm;0  WITH  A PARTICULAR  ACLS  ORIENTATION 

C a're  CALCULATEO 

C CALCULATE  TRAMSC»NP'MT4LS  ONLY  ONCE 

CSSCS*COS(»MlEl*3lN(TMrT4EI»SINI<I£>-COS<SI£>»SIh(PMlC» 
CPCT»COS<FhIE)»CCS(ThETA£»  


C clear  total  FORCES  AMO  TORQUES  TO  EERO 

fORCT*0.0 

lTPX«b,0 

TTPZ»0.0 

TCPX*b.O 

TCPI<0.0 • 

TORF7.0.0 

TOROTX.O.O  ' 

TORaT7»0,  0 SPORCEtVo.  JTOFJUE'xpO.  STCRCOEZ»0. 

C FORCES  and  TOTOUES  INPEPENOtNl  OF  JtGM£NIS_lNCIVlOUAlLTL 
C HEAVc,  FORCES  CUSHION  AND  TRUNK 

FCP»PCH»ACM  _ _ _ ; 

FTPaPTK'ATKCN 

jl_ compute  velocity  for  drag  forct 

V*V£L«*CSS;S*SINKRT*CPCI 

_ SICN«1,  ^ _ 

IKV.GI.0i3i  SICN.-l, 

c meavl  drag  force  _ 

FOF«  0.  S*HOC*»UA‘RhO*  V*V»SIGN 7 

_CJ)RAC  TOROUE  

iof7*fcf»c-:mfx  ■ 

TOFxx-roF •CFNFZ 

'C  FOPCtS  ANO  TORIUES  OCPENDTHT  ON  SEGNCMTS  INDIVIDUALLY 
C SUM  IROIVICIOUAL  SIGmpnTS  TO  FIND  TOTALS 
CD  101  IK.NSTOP 
C CUSHION  PRESSURE  TORQUES 


Rra  A.  -TIIABIE  COPV 
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TC*>r«tCP?*  ixc**e  uchui  i-*cmrm  1 1 

T1PJ»1T“2*  «4T<C».IIl  (•/ItKCHSCIU  ) 

TTP««1T»X- II 
in  ittCN!  It  I.  GT.  0.  i.cR.  («rKC’iR(ii.‘;T.  0. 1 1 <io  to  iii 

CO  TO  101 _ 

lu"  VCLTi  I rN<*T»CPCT»0»Ml»l*ti(lll-CCI*OTHtT*»(iTk(;|.Pri 

f CPC»-V£LT»0*P»J»P£RHI  I _ _ 

fORCT.»ORCT*rORO  . - . 

TO<tOTl.TaPlTl*(«TKIII-CCl*PORO 

TOROTXiTOP.lTX- <1T<(  Ii.f<-1  •roRO 

invClx.EO.O.OI  CO_TO_101  _ 

TORP?«TORf »-ITC3*CPCT‘  <PT<* lirkcNI I II «AT<CNRI 1 1 1 ‘U 

tot  CONTINU-  _ 

C COMXiTlON  CP  P»'c€  ANfl  fORQue  COIIPOMeHTS 

C total  heave  potcs  

fORC£x«(PC»*PT»*PORCT«rori»CPCT 

. _ rT*AHA«imo.CCTtO.  I ■ 

C total  TOaouE  X AXIS 

TOROy£X«TC»X«TTOX»TOP.OTJ5»Tp^* 

TI'TOROUEX 

CJOTAL  TOOOUE  1 AXIS  

TOROUIl«TC’2*TTPl»fOROT2»TOr2»T6»i>T 

T2»T0RCyi2 : 

C STATc.  EQUATIONS 

C»»*THi  STATE  VAIUSLES***  

C llP»LH..PLfxUN  PRESSURE  (CACEI 

C 2IPCH.. CUSHION  PRESSURE  (CACEI ^ ; 

e l|PTX..TaUN<  PR-SSURE  (CACEI 

c nsiNXRT,. vertical  sink  rat:,  positive  upraros 

C 5ITCC..CC  e(.EVAT10N 

C 6I0»MI,  .PITCH  RATE. VEHICLE  PRA  l£ ; 

e TIOTHfTA. .ROLL  RATE.  VEHICLE  RRAHE 

C »|Th£T4£..£U.XRI1!J  roll  AN'.L- 

c oiPH:E..tuLiRHH  PITCH  angle 

C_1tlSl!..EULERIW  TAH  ancle  (APPROX. _;£ROJ ; 

e lltXV.,5ISPLN  or  »r:ssu®e  RSlUP  valve 

C t2ivv..v:ic:iTY  or  PRESSURE  RELltr  VALVE 

C 13IOrtNx..rAM  AIR  INERTAKTE  riOH 

OCRT(tl«(CX<*(“RLH.-»AT|/VPLH|»(OrAHX-QRLTKI  

CfRTdl  • (C<K«  (PIK  ♦PATl/VTXI  •(OP11X-OIXCH-OTKAT*OVTKI 

C. CUSHION  PLOX  AROVE  C«OU'(D  -pFECT  TRANSITION  20NE 

OCNf  IpOPLCT'OTXCm-OVCm  ’ 

e CALCULATE  CROUXO  E»*ECT  TPAKSIJION  19^1. 

T90UNO«OC‘HY1MI.  .SEC! 

_ ieouNO^cC'Hri  

C OETECNINE  IP  a;ls  in  transition  ione 

ir(TCO.CT.TTO'JNOICO  Tp  13 ^ 

ir (rcc.cT.BOouNoi  CO  to  in 

CO  T0_16  ^ 


R*OVt  TRANSItnN  70ii 

t3  OCNATiOCNrT 

IFLAC.e 

lOPO 

_ N0«0  » 

CO  TO  15 

IM  transition  ione 

_1C_  IFLAC.I 

NOalOO 

10«ITIX (AOS  ((TS3UNP-»CCI/(I10UNO-nROUNOI*rLOAT (NOI  I I 

irdO.CE.NSlIO'NO 

CO  TO  12  _ 

c 

C In  CNOUNC  EFRECT  70nE 


V 


95 


AFIT/GE/EE/77D-43 


r 

t 

i! 


I 

I 


I 


Ik  NOd 

IO«l  . 

jrL*C»? 

C COMPUTE  CU'MiOi  TO  ITMO^phTRE  ELOW  

ir  OCMiTtPLOl  T(HO-tOt  /'LO»T(NOI*OCHET»piOAT  t lOI  /float  (NO)  •QCHkT 

C CUSHION  PPiSSUT-:  DERIVATIVE  

*S  OC«V  (?)*  IOPLCM«'JTVCM.cCHAT-avCH»CvCHF»PCM»OEA»l  J)/(PTK»PTICI  )/ 

I (VCM/ICV<* l“CH‘®A tl I »PVShP/PTiO  

C pUSMlCN  pscSSUR:  is  zero  AOOVE  transition  ZONE 

iniFLAO.ET.O.OI  PCH-.O.fl 

kk  RETURN 

£NO ^ 

PARAMETER  VAluES>IOIAV>0 

_ OMAXiZO. 

MSS»TT.5  CC«-l.lT 

CO*  l.S  ,»F«0.  ,OHA<l$S.  ,jiOC«l,  ,L5-<i.ir5  tO»>Al7  -A«.&RS 

£«.SZ<>.L«S.  IT  .HVIsL.OO  I AM>3.NE*L  .SHi.083I  .LOO.AO 

Air..JZ2  _ 

APC>0.  .APT'. 695  .APA>0.  .AAT>.796  ,VCD*0. 

DPV«e^ FAT.2I16.*  ,T£M*70.  _ 

.CVi(<l.i.  ,CPA»C.  .CAP^.S  .CPC*0.  .CPTs.S  .CTC«.9 

_ CTA*.<.  ,CCP«.5  ,0£C*.l  CP»«  1.17  ,Cr2«0,  .CPC'S. Z 

*A  LO'O.  .XU  tO'l.  .ZA  LP»0.  .HAttO'O.  • .MU  L0»0.  ,XB)L6»-'iz 

Y90LD»Q.  «xR  LO'Sji lTO»t-0  = 0, tS-ROl-PiO. 

10  LO'l'.  '.NeSLO.n,  .NS  LC'O.  .NOALO'b.  .RUOLO'l.  .8  10*19. 

TCOEN*a.<.  .THREN.ZI53.  .XO  FN*-b.  1 . GA»EN*1.  _ .GAZEN*0. ,P0  TF.155^^ 

PI  TF*l/.‘.6.Z0  TF,o.  ,Z1  ’F«t55.7S 

.Z«  TF  ZMSa.^  

Zl  TF  2.0..P0"  TF  2*151.  TS.pT  t>  2*iT.v6, 

6K1IT  S»l.  _ .GNLIT  S*1. 

'ci  MC  IPI.Z.C'Z  M;  HJ.75.C3  HI  1*2. 5. G*.  MC  1»0..C>(IIT  •.•«1..6<LIT  7*1. 
CKIIT  A*-I..CKLIT  8*1. .Cl  "C  Z*1.JA.C2  MC  2«5.7S.C3  MC  Z*2.5.CA  MC  Z«0.. 
CRlIT  5>l..CKLtT  SH.  .GRIITia**!.  .GklITIO*!.  .ALSAV'O.  .S  AV*76. 

UH  AV«0.  .WV  AV'.O, ,MM  AV*0,  _.PH  AV*0.  .OMIAV.O.  .RNiAV*0.  

MAitO»77.9.C  LCtv.  ,G<LIT  t«l.  .GKIIT  1*1.  . IXXSC* ZNiv, , I TYSO* 1 BfcO.  . 

IZZ5O«.oeZ.,IXZS3*.2oa._,T0  en*q.  ,to  to«o. 

FTtS3>0.  .»Y2S1»0.  .TZZS3.J,  .fx2S3*0,  .TXlSS»().  ,TZ1SS»0.  , 

FX3$3<0.  .FriSI'T.  .FZ333.0. 

XP  L0.-99.6?Z.rA0L9*-  .123ft6.YP  lO*-. 035592720  LO*-5,P257, 

MAOLO*-. 75599. YR  LD*.  12081  .1®  L 0»- 1 1 3. 7i. XOILO*-.  055  0 52. MO  L0*-l75.86j 

YORLO*350,  17,lR  10»- I*.  3 15.U»  ID*. 1 1. 8* 1 . ZP  10*- 1 0. 957. ZCELO*- 1 5.  5 52. 

IORLO*607,  .NS  10*100. 93. HO  10*13.  72 5 . H0£10*-86S. 1 3,L0Air»-3. 579 5 . 

'n0R10>-119. 

INITIAL  CONOITIO'IS*  __  

PPLFM.137,91.PTi<*M*125.  05.PChFM*0,  ,Th  EN>1515.8.U  SD>17.  033, 

V S0*6.S2S5.h  <0*0. ,P_  S9*0. lO  SO*0  .0 R Sl)*0.j 

R0LS0*0.  .PirS9*-l.P  .YAMSd'*0.  .ALISD*2.5  0.XI~TF«0', 

X2  TF*0.  .X2  IT  7*0.  .12_3’  »'0.  _ _ -X?  12  I.P, 

X2  IT  1*0. .xf  TF  2.0,.X2~TF  2*0. ,X2  ir'9*6..X2  ITIO'O. 

_P«IINT  CONrP31«..T_IN2*0.  ^ .THAX*10,jOuTRAII*1 . PRATE. _l,SIMUiU£ 


I 


i 

1 


AFIT/GE/EE/77D-43 


Vita 

Edward  Arthur  Kenney  was  bom  18  July  1948  In  London,  Ontario, 
Canada  to  James  A.  and  Evelyn  V.  Kenney.  After  graduating  from  G.  A. 
Wheable  Secondary  School  In  1967,  he  entered  the  Royal  Military  College 
(RMC)  In  Kingston,  Ontario.  He  graduated  from  RMC  In  1971  with  the 
degree  of  Bachelor  of  Science  In  Electrical  Engineering,  and  a commis- 
sion In  Che  Canadian  Armed  Forces.  After  graduating  from  the  Primary 
Flying  School  and  attending  the  Flying  Training  School  at  CFB  Moose- 
Jaw,  Saskatchewan,  he  reclassified  to  be  an  Aerospace  Engineer.  Sub- 
sequent tours  of  duty  included  CFB  Shearwater,  Nova  Scotia,  as  the 
Deputy  Maintenance  Records  Officer  and  Heavy  Maintenance  Repair 
Officer;  H.M.C.S.  Preserver  as  the  Fleet  Air  Maintenance  Officer; 
and  National  Defence  Headquarters,  Ottawa,  in  the  Directorate  of 
Avionics  and  Armament  Subsystems  Engineering.  Cape  Kenney  was  as- 
signed to  Che  United  States  Air  Force  Institute  of  Technology  in 
June  1976  in  the  Graduate  Guidance  and  ControJ  Program, 

Mailing  Address:  National  Defence  Headquarters 
Ottawa,  Ontario  Canada 
KIA  0K2 

Attn:  DPCO/AERE 


97 


AFIT/CE/CE/770-43  in:CL\SSIFIED 


• eCOWTV  CLASSiriCATlON  OF  THIS  PAGE  flWiAO  Dal*  Enimtmd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


ue. 


IT/GE/EE/77D-A3 


I*.  COVT  accession  no 


J.  RECIPIENT'S  CATALOG  HUMBER 


TiTLS  fmtrsctnri 


[6  Design  and  Simulation  of  a Takeoff 
Stabilization  System  for  an  Aircraft 
ith  an  Air  Cushion  Landing  System  • 


n 


a 


I Tf-ti  wr  nbf'vni  • r'&nivw  W’ 


t.  performing  org.  report  number 


7.  AUTHORfa; 


•.  contract  OR  grant  NUMBERr*; 


Edward  A. j Kenney 


9.  performing  organization  name  ano  address 


Air  Force  Institute  of  Technology  (AFIT/EN) 
Wright-rPatterson  AFB,  Ohio  A5433 


10.  PROGRAM  element,  project.  TASK 
AREA  B work  unit  NUMBERS 


I).  controlling  office  name  aho  aooress 

Vehicle  Equipment  Division,  Mechanical  Branch, 
(FEMB) , Air  Force  Flight  Dynamics  Laboratory, 
Wright-Patterson  AFB,  Ohio  A5A33 


U.  monitoring  AGENCy  name  a AOORESSri/  <D//*r*nl  /rom  ControlUng  OtIIca) 


IS.  securityc'l'aSH’oI^IW/Ifo^ 

Unclassified 


tSa.  DECLASSIFICATION/OOWNORAOINC 

schedule 


U.  DISTRIBUTION  STATEMENT  (ol  ihim  Rapari) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  kbotroef  ontorod  In  Block  30,  II  dificrcnr  Irom  Rmport) 


ic  release;  lAW  AFR190-17 


JER^L  F.  wtoSs,  captain,  USAF 
Director  of  Information 


19.  KEY  WORDS  (Continuo  on  rororoo  otdo  II  nocoocmry  md  Idontlly  by  block  number) 

ACLS  ^SY 

SAS  acts 

Coanda  Effect 
Thrust  Vector  Device 
Aircraft 


m 


ABSTRACT  (Continue  on  reeeree  elde  II  noceeeory  m*d  Idontlly  by  block  number) 
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J As  well , the  thrusters  would  be  set  forward  of  the  aircraft  cen- 
ter of  gravity  and  could  be  activated  in  tandum  to  aid  in  pitch 
control . 

The  Jindivik  Remotely  Piloted  Vehicle,  an  Australian  target 
drone,  vras  fitted  with  an  ACLS  and  taxi  tests  showed  the.  insta- 
bility and  need  for  a stabilization  system.  Subsequent  use  of 
Jindivik  wind  tunnel  and  taxi  test  data  served  as  the  basis  for 
the  development  of  the  roll/pitch  control  system  presented  in 
this  paper.  Due  to  computational  problems  with  the  air  cushion 
model  of  the  computer  program,  the  controller  designs  could  not 
be  completely  verified;  but  expected  trends  in  pitch,  roll,  and 
yaw  control  were  shown. 


